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ABSTRACT 

We investigate how "extra" or "excess" central light in the surface brightness profiles of cusp or power-law 
elliptical galaxies relates to the profiles of ellipticals with cores. The envelopes of cusp ellipticals are estab- 
lished by violent relaxation in mergers acting on stars present in gas-rich progenitor disks, while their centers 
are structured by the relics of dissipational, compact starbursts. Ellipticals with cores are formed by the subse- 
quent merging of the now gas-poor cusp ellipticals, with the fossil starburst-components combining to preserve 
a dense, compact component in these galaxies as well (although mixing of stars smooths the transition from 
the outer to inner components in the profiles). By comparing extensive hydrodynamical simulations to ob- 
served profiles spanning a broad mass range, we show how to observationally isolate and characterize the relic 
starburst component in core ellipticals. Our method recovers the younger starburst population, demonstrating 
that these dense concentrations survive spheroid-spheroid mergers and reflect the degree of dissipation in the 
initial mergers that formed the penultimate galaxy progenitors. The degree of dissipation in the mergers that 
produced the cusp ellipticals is a strong function of stellar mass, roughly tracing the observed gas fractions 
of disks of the same mass over redshifts z ^ 0-2. The strength of this component strongly correlates with 
effective radius at fixed mass: systems with more dissipation are more compact, sufficient to explain the dis- 
crepancy in the maximum phase-space densities of ellipticals and their progenitor spirals. The survival of this 
component together with scattering of stars into the envelope in re-mergers naturally explain the high-Sersic 
index profile shapes characteristic of very massive core ellipticals. This is also closely related to the kinematics 
and isophotal shapes: only systems with matched starburst components from their profile fits also reproduce 
the observed kinematics of boxy/core ellipticals. The final "core-scouring" phase of core formation occurs 
when a black-hole binary formed in the merger scatters stars out of the innermost regions of the extra-light 
component. It is therefore critical to adopt a physically-motivated profile-decomposition that accounts for the 
fossil starburst component when attempting to quantify scouring. We show that estimates of the scoured mass 
that employ single-component forms fitted to the entire galaxy profile can be strongly biased. 
Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation — 
galaxies: nuclei — galaxies: structure — cosmology: theory 



1. INTRODUCTION 

Early work to incorporate gas physics and star forma- 
tion in models of the formation of bulg es and elliptical 
galaxie s through mergers demonstrated (iMihos & Hernquist 
Il994allbll3) that dissipation in disk-disk mergers leads 
to central starbursts lik e those observed in UL IRGs 
(lJoseph& Wright! Il985h . iMihos & Hernquist! (! 1994a) pre- 
dicted that this process should leave an observable signa- 
ture by imprinting a two-component structure into the surface 
brightness profiles of merger remnants. That is, the com- 
pact stellar remnant of the starburst should have a steeper 
density profile than the inward extrapolation of the outer 
Ide VaucouleursI (! 19481) r'/'^-law shape of the main body of 
the galax y. The "extra light" in the starburst remnant is 
required (!Hernquist et alJ !l993!) to explain the high mass 
and phase space densities of bulges and elliptical galax- 



ies ('Ostriker'1980'; 'Carlberg T986^ ^unh 1987; 'Kormendy 
1989) and their fundamental plane correlations (Dressier et al. 
1987; Diorgovski & Davis 1987; Kormendv &Sandersll 19921; 
Rothberg & Joseph,2006- Hopkins et al...2008blK 
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When predicted by lMihos & Hernquisl (!l994a|), ext r a ligh t 
was not known to exist in normal ellipticals. !Sersid (!l968!) 
functions were replacing r'/"^ laws as standard machinery to 
interpret brightness profiles dCaon et al.l |1993 |). The steep- 
est central pr ofiles were calle d "power laws" (jperra rese et al] 
1994; Kor mendy et al.1 !l994t !Lauer et alJ !l995!: !Byun etaP 
,1996.) . Both descriptions could be taken to imply, absent cores 
or nuclear star clusters,^ that the brightness profiles of ellipti- 
cal galaxies are almost power laws with no breaks indicative 
of two-component structure. However, Kormendv ( 1999) dis- 
covered in a small sa mple of normal ellipticals - and later 
!Kormendy et al.! (!2008i) . in all the known low-luminosity el- 
lipticals in the Virgo cluster - that the inner "power-law" 
or "cusp" indeed represented a central "extra light" compo- 
nent - the inner profile of a dense, centrally concentrated 
rise above the inwa r d extr apola tion of an outer S ersic pro- 
file. iFerrarese et"a!] (l2006h and ICdte et al.1 (|2007|) identified 
qualitatively similar features - "luminosity excess" in the cen- 

^ It is important to distinguish nuclear clusters (stellar "nu clei") from ex- 
tra ligh t components - we demonstrate the clear distinction in lHopkins et all 
I2008al) . Nuclei are extremely small (typical radii ~ 1 - lOpc, and mass 
fractions ~ lO^^Mg^i) relative to starburst/extra light components (radii ~ 
0.2-2kpc and mass fractions ~ 0.1 Mgai). The structural scahngs and funda- 
mental plane relations of stellar nuclei are similar to those of globular clus- 
ters, nearly perpendicular to and with an opposite physical sense from extra 
light components. Likewise their stellar population ages, metallicities, and 
kinematics are very different from extra light and elliptical galaxies I Caroll3 
1 1999; Boker et aL2004;,Walcher et aL200fc.C6te et al..200fc,Hopkins et all 

boosj) . 
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tral regions of ellipticals (for a detailed comparison of these 
features, see the discussion in iHopkins et alJl2008ah . Exten- 
sion of this to the profiles of recent gas-rich merger rem- 
nants observed by iRoth berg & JosephI (2004) showed simi- 
lar signatures. The "power law" profiles identified by pre- 
vious authors form a part of this extra light; the conclusion 
that there are two components, however, did not become 
clear until observations from a variety of telescopes made 
it possible to create profiles with sufficiently large dynamic 
range. These authors suggested that the observed extra light 
is the signature of merger-in duced starbursts, as predicted by 
iMihos & Hernquisa (Il994ah . and it is now becoming empir- 
ically established that essentially all "power-law" or "cusp" 
ellip ticals show thi s behavi or 

In iHopkins et al.l(l2008dl) (hereafter Paper I) we showed that 
observed extra light can indeed be identified with the central 
density excess produced in simulations of gas-rich mergers. 
We developed a formalism to fit surface brightness profiles 
with two-component models in a manner that accurately re- 
covers the physical decomposition between a central dissi- 
pational component (the remnant of the merger-induced star- 
burst) and an outer dissipationless component (the product of 
violent relaxation acting on stars formed in the pre-merger 
disks before the coal escence of the mergin g galaxies). We 
used this machinery in lHopkins et al.l (l2008ah (hereafter Paper 
II) to analyze observed "cusp ellipticals"^. We showed that 
this method could be used to recover the dissipational com- 
ponents in cusp ellipticals, and that such components were 
ubiquitous in the local cusp ellipti cal population (in sam- 
ples of ~ 10 0s of such galaxies from lLauer et al.l ilOOl d) and 
iKormendv' et al. (2008)). In essentially every case, simulated 
dissipative merger remnants provided good matches to the ob- 
served profiles of cusp ellipticals. Mutually consistent decom- 
positions of the observed galaxy profiles and the simulated 
merger remnants demonstrated that the structure of cusp el- 
lipticals as a function of mass corresponds to that predicted 
for the remnants of gas-rich, spiral-spiral mergers with real- 
istic properties. That is, the degree of dissipation needed for 
mergers to explain the densities and scaling laws of ellipti- 
cals corresponds well with our empirically recovered starburst 
components and agrees with the gas fractions available in ob- 
served progenitor disk galaxies of appropriate masses. 

However, it has been argued that the situation may be 
different for more massive ellipticals, which appear to ex- 
hibit central "cores," and are known to display boxy isopho- 
tal shapes and slow rotation in contrast to the disky isopho- 
tal shape s and rapid rotation seen in le s s massive cusp el- 
lipticals (iKormendv & iningworthlll982t IPavies et all 1 1983 



Davis et al] 



Bender et al 



19851 Mrzeiewski et al. 1987; Bender 1988 
119891 ri992; ' Peletier et al.nl990i) . Cores were 



first seen in ground-based observations of nearby, high- 
luminosity ellipticals as central regi ons of nearl y constant sur- 
face brightness (.King. 1978; .Young^etaLlfTgTi lLauer.l985bl; 
lKormendvlll985ail987ah . They were cuspier than isother- 
mal cores, but the functional form of the density profile as 
r ^ was unknown. Later, HST images showed that nearly 
all ellipticals have singular starlight distributions in th e sense 
that the surfac e b rightness diverges a s S](r) ~ r~'^ (Lauer et alJ 
[19911 fT992allbt ICrane et al.1 fl99l IKormendv et all 119941: 

* The term "cusp elliptical" is used in various ways in the literature. Unless 
otherwise stated, we use this name to refer to ellipticals that do not have 
a central resolved core. In pubHshed papers, these are called "power-law 
ellipticals", "coreless ellipticals", or "extra light ellipticals." 



iFerrarese et al.l [19941 ; iLauer et alj|1995h . In low-luminosity, 
early-type galaxies, 7 typically decreases only slowly as the 
center is approached, and a steep 7 > .5 cusp continues in to 
the HST resolution limit. ^ Lauer et al] ([1995 ) classified these 
as "power-law" galaxies. In more luminous ellipticals, the 
steep outer density profile shows a r obust break to a sh allow, 
inner power-law cusp with 7 < 0.3 dLauer et al.lll995h . The 
"break radius" r/, corresponds to th e "core radius" T c mea- 
sured in ground-b ased observations (IKormendv et al.l fr994). 
iLauer et al] ( 119951) continued to call these "core galaxies," 
even though the shallow cusps in projected brig htness imply 
steep and singular cusps i n luminosity density (ILauer et alj 
fT995HGebhardt et al.l[T996l) . 

The division of central structure into two families was 
motivated by the observed bimodal distribution of cusp slopes 
7 (Gebhardt et al. 1996; Lauer et al. 2007a), but this aspect of 
central slopes participates in a larger, longer-recognized divi- 
sion of the elliptical population: the typical giant, core ellipti- 
cal has different physical prope rties from the typical normal 



luminosity, cusp elliptical (e.g. Davies et al. 1983; Davis et al 
1985; Jedrzeiewski et al. 1987; Bender 1988; Bender etal 
1989, 1992; Peletier et al. 1990; Kormendv & Bendei 
1996; Fab eretal.. .1997.: .Simie n & Prugn iel ',1997a.biii 
Emsellem et all l2004f 120071: [McDermi d et alj l200i 
Cappellari et alj 



Massive giant ellipticals rotate 
slowly and have "boxy" (rectangularly-distorted) isophotal 
shapes, characteristic of systems supported by anisotropic 
velocity dispersions (ISchwarzschildl Il979t IdeZeeuwl [1981 
lBinney&Tremaindll987h : iFaberetalJ (Il997l) showed that 
these relate to the observed "core" population. In contrast, 
less massive ellipticals (and SO galaxies), where the power- 
law population predominates, rotate more rapidly; they have 
more isotropic velocity dispersions and look like they have 
embedded disks ("disky" isophotal shapes). 

These differenc es thus led naturally to the idea, de- 
veloped,^n e.g. lFaberetal.1 ( [1997): Kormendv (1999); 
Guillen etal. (2000); Rest et al. (2001); Ravindranath et al] 
(2001) : jLaineetaL (2003,); .Lauer etal.i (2005, 2 007al): 
Fen-ar ese et al.l (l2006l)nCdte et alj (l2007l) : [Kormend v et alj 
(i2008t and references therein), that disky, rapidly rotating 
cusp ellipticals are direct products of gas-rich ("wet") merg- 
ers, whereas boxy, slowly rotating core ellipticals have been 
shaped by subsequent dissipationless ("dry") re-mergers of 
two or more (initially cuspy gas-rich merger remnant) ellip- 
ticals. Several questions therefore arise. How did core ellip- 
ticals form? What were their progenitors? It has been shown 
that mergers of bulgeless disks fail to reproduce the shapes 
and k i nematic properti e s of these galaxies (see e .g. iBarnes 



ana Ki nematic properti es or t nese galaxies (see e.g. marne^ 
[19881 lHernauisa[l99l [1991 iNaab et alJ [2006i ICox et al 
I2OO6F). Furthermore, if disks are their progenitors, then these 



" progenitors, tnen these 
systems would not be able to avoid dissipation, because spi- 
ral galaxies contain gas. If lower-mass cusp ellipticals are 
the progenitors of core ellipticals, what happened to the extra 
light in those cusp ellipticals? 

Nearly all numerical experiments find that light profile 
shapes ar e, to lowest order, preserved in dissipationless merg- 
ers (e.g., Bovlan-Kolchin et al.l 1200 5). So, core ellipticals 
should be "extra light" elUpticals in the same sense as cusp el- 
lipticals - i.e. their profile reflects a combination of an outer, 
low phase-space density violently relaxed remnant of stellar 
disks, and an inner, compact component originally formed via 
dissipational processes (in whatever process formed the pro- 
genitors, that would allow gas to lose angular momentum and 
reach these densities in the first place). If their last merger 
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was dissipationless, is the amount of dissipation in the "orig- 
inal" spheroid forming merger (i.e. that which formed their 
progenitors, if these are re-merger remnants) important or rel- 
evant to the properties of the z = galaxy? Is any memory of 
the merger history preserved, and is there any way to obser- 
vationally recover this information? 

There have been suggestions of interesting behavior in the 
observations, but they have largely lacked an interpretive con- 
text. Models of core galaxies generally presume that the pro- 
jected stellar brightness decreases monotonically outside the 
core with a monotonically increasing logarithmic slope - the 
fact that many core ellipticals can be reasonably well-fitted by 
such descriptions has made the multi-component structure of 
such objects more ambiguous. A closer look at the existing 
data, however, shows evidence in some galaxies of charac- 
teristic features such as inflection points in the slope outside 
the core, similar to a "smoothed" version of the features of- 
ten seen in cusp ellipticals where the profile transitions from 
being dominated by a dissipational ("extra light") component 
superimposed on a background dissipationless outer envelope 
(see e. g. RatcHff & Burstei n 1982; Barbon et al. 1984; Laugi 
Il985bl: iLauer et al.l l2007al; [Kormendv et all l2008h . More- 
over, even where profiles are smooth and monotonic, the in- 
ner portions of core galaxies typically rise well above an 
r'/'* law fitted to the envelope. Instead, single profiles fit- 
ted to the galaxies are forced to higher Sersic indices that 
rise more steeply at small radii, reflecting and including 
the dense, high surface brightness central extra light com- 
ponent, and yielding the high Sersic indices (rts ^ 6 - 8) 
characterist ic of core elhpticals fitted in th i s manner (e.g. 
Ift-uaniel & Simienlll997t iTruiillo et al.ll2002y Ferrarese et al.l 
|2006; Kormendv et al. 2008). Cote et al. ( 20Q3) point out that 
there appears to be a smooth transition from low-nj outer pro- 
files with a steep central rise above the corresponding inward 
extrapolation, to high-n, outer profiles with a corresponding 
rise implicit in the Sersic profile, and hence not explicitly ap- 
pearing above this threshold. Without a physical motivation 
for decomposing this inner rise and outer envelope, interpre- 
tation of this phenomenon has been restricted to the empirical 
notation of the best-fit Sersic indices. In the context of the 
present work, however, it highlights the potentially composite 
nature of core galaxies. 

It is generally believed that the connection between the 
merger history of galaxies and their nuclear profile slope 
("cusp" or "core," on scales much less than the effective radius 
or the scales of the extra light) arises because of "scouring" 
by a b in ary black hole (for a rev iew, see lGualandris & MerrittI 
l2007h . iBegelman et al.l (Il980l) first pointed out that binary 
black holes coalescing in a dissipationless galaxy merger stall 
(i.e. are no longer efficiently transported to the center via dy- 
namic friction) at radii '--^pc, larger than the radii at which 
gravitational radiation can efficiently dissipate energy and 
merge the binary - the so-called "last parsec problem." They 
noted that significant gas content can provide a continuous 
source of drag and friction and solve this problem in gas- 
rich mergers, but that in "dry" mergers, the binary will remain 
stalled for some time, and will harden by scattering stars in the 
nucleus in three-body interactions. This will continue, flatten- 
ing the nuclear slope, until sufficient mass in stars Mbh, 
by simple scaling arguments) is ejected to merge the binary. 
It is therefore of particular interest to estimate the stellar mass 
which must be scattered to explain the slopes of cores, as a 
test of scouring models and (in such models) a probe of the 



galaxy merger history. However, such estim ates have been 
ambig uous (and often co n ti-oversial: see e.g. iFerrarese et al. 
2006; Laueretal.ll2007allbl; ICote et al.ll2007t iKormendv et al. 
2008), in large part because of the lack of an a priori physical 
model for the profile shape. Understanding the global profile 
shapes and extra light in core ellipticals is therefore impor- 
tant to reliably estimating how their nuclear profiles have, in 
detail, been modified by scouring. 

It is also important to recognize that there can be both 
continuity and bimodality in the cusp/core populations. Be- 
cause the expected number of major mergers in the forma- 
tion history of a typical elliptical is not large (~ a couple), 
it should be a relatively Poisson process: a significant num- 
ber of ellipticals (especially at low masses) will have expe- 
rienced only the original, single major gas-rich merger that 
transformed them int o ellipticals since z '- ^ 2 - 3 (see e.g. 
Mailer et al. 2006; Hopkins et al.l I2008ell3 ISomerville etaU 
i2008i;iLin et al.ii2008i) ; others (especially at high masses) will 
have experienced ~one or two subsequent major mergers, 
which will tend to be "dry" spheroid-spheroid mergers. Al- 
though there might be some intermediate cases, to the ex- 
tent that properties (such as kinematics and isophotal shapes) 
are affected by the last major merger, there should be signif- 
icant differences between those whose last merger was dissi- 
pational (with a disk that contains some mass in gas) or dissi- 
pationless (spheroid-spheroid). 

However, although the last merger may be dissipational or 
dissipationless, the total amount of dissipation in the forma- 
tion history - the mass fraction formed dissipationally, from 
gas losing angular momentum in mergers and participating in 
nuclear starbursts - should be continuous across either pop- 
ulation (for e.g. a given mass and original formation time). 
Although some objects may have had subsequent dry merg- 
ers, they still formed stars in a dense central concentration 
in the original, gas-rich merger that formed the progenitor 
spheroid; this process will be the same regardless of whether 
or not the system is destined for a future dry merger. As 
a function of mass, this should broadly reflect the gas frac- 
tions of ultimate progenitor disks at the spheroid formation 
times, and as such is a continuous function of mass, star for- 
mation history, and formation time. This continuity in dis- 
sipational content, to the extent that it effects the structure, 
fundamental plane correlations, and stellar populations of 
spheroids, is reflected in the continuity of e-g- t h e fundamen- 
tal plane (recently, s ee 'Cappellari et al.J 120061; |Bolton et alj 
120081; iHvde & Berna rdi 2008), the stellar a ge and m etallicity 
versus mass rela tions (Traaeretal. 2000; N elan et al.l 120051: 
iThomas et al.l2005: Gallazzi et al. 2005, 200 6}), and the color- 
magnitude relation dStratevaet al. 2001; Bal drv et al] 12004). 
As there has been considerable observational debate regard- 
ing the degree of continuit y or bimodality betwe en cusp and 
core populations (see e .g. Ferrarese et al. 2006; Laue r et akl 
I2007al; ICdte et all 120071: Eormendy et al. 2008), it is clearly 
of interest to identify properties that are or are not expected to 
be continuous across detailed merger and re-merger histories 
within the spheroid population. 

In order to understand the structure and formation history of 
the core elliptical population, we therefore extend our study of 
merger remnants and cusp ellipticals in Paper I and Paper II 
to the core population in this paper. We wish to test the hy- 
potheses that these systems were, in fact, originally formed 
(i.e. their progenitors were formed) in gas-rich mergers (al- 
beit potentially modified in gas-poor re-mergers), and that the 
original degree of dissipation can be empirically recovered. 



4 



Hopkins et al. 



and is the critical parameter that can explain their densities, 
scaling relations, profile shapes, and sizes. 

In § 12] and § [3] we describe our set of merger simulations 
and the observational data sets we consider, respectively. In 
§ m we study how light profiles of gas-rich mergers, which 
we studied in detail in Paper I and Paper II, evolve in sub- 
sequent re-mergers of such ellipticals, and demonstrate that 
our fitting procedures designed to recover the original dis- 
sipational/starburst component in gas-rich merger remnants 
can be applied to re-merger remnants. In §|5]we investigate 
how gradients in stellar populations, imprinted by the extra 
light and original gas-rich merger, are affected by re-mergers. 
Readers interested primarily in our comparison of the prop- 
erties and scalings of dissipational components in observed 
core ellipticals may wish to skip to § |6] where we compare 
our simulations with and apply our fitted galaxy decomposi- 
tion to a wide range of observed systems. In §|2]we use these 
comparisons to study how structural parameters of the outer 
stellar light and inner extra light component scale with galaxy 
properties, and compare them with the extra light components 
in gas-rich merger remnants and cusp ellipticals, and examine 
how the existence and strength of the extra light component is 
related to galaxy structure and drives galaxies along the fun- 
damental plane. In §|8]we consider the global isophotal shapes 
and kinematic properties of re-merger remnants and how they 
depend on extra light content. In § |9] we examine how re- 
mergers and issues of profile fitting relate to the outer Sersic 
profiles of core ellipticals and re-merger remnants, and com- 
pare results obtain with different choices of empirical fitting 
functions. In §[TO]we demonstrate that this can affect esti- 
mates of "missing light" on small scales, and explain how this 
effect arises and what it means for a proper physical under- 
standing of nuclear light profiles. Finally, in §[TT]we discuss 
our results and outline future explorations of these correla- 
tions. 

Throughout, we adopt a JIm = 0.3, ft\ = 0.7, i/o = 
70km s"' Mpc"' cosmology, and appropriately normalize all 
observations and models shown, but note t hat this has little 
effect on our conclusions. We also adopt a IChabrierl (|2003|) 
initial mass function (IMF), and convert all stellar masses and 
mass-to-light ratios to this choice. The exact IMF system- 
atically shifts the normalization of stellar masses herein, but 
does not substantially change our comparisons. All magni- 
tudes are in the Vega system, unless otherwise specified. 

2. THE SIMULATIONS 

Our merger si mulations are described in Paper I and Paper 
II rfor details. seelSpringefcOOSU Spri ngel & Hernquis3l2002l 
I200I ISpringel et al.T]2005[r They are fully hydrodynamic 
simulations of galaxy-galaxy mergers, including physical and 
empirically motivated models for black hole accretion, star 
formation, a multi-phase interstellar medium, and feedback 
from supernovae, stellar winds, and black hole growth (we 
show in Paper II that our conclusions are robust to variations 
in these prescriptions). The typical gravitational softening is 
^ 20-50pc O.QlRe), and hydrodynamic gas smoothing 
lengths in the starbursts are smaller, which we show in Pa- 
per I and Paper II is sufficient to properly resolve not only 
the mass fractions but also the spatial extent of the extra light 
components of interest here. 

We consider a series of sev eral hundred simulations of col- 
liding galaxies, d escribed in iRobertson et al.l (l2P06a b) and 
ICox et al.l(l2006abt) . We vary the numerical resolution, the or- 
bit of the encounter (disk inclinations, pericenter separation). 



the masses and structural properties of the merging galaxies, 
initial gas fractions, halo concentrations, the parameters de- 
scribing star formation and feedback from supernovae and 
black hole growth, the equation of state of the interstellar 
medium gas, initial merger redshifts, halo virial velocities, 
merger mass ratios^, and initial black hole masses. 

To this, we add a subset of spheroid-spheroid "re-mergers," 
representative of gas-poor or "dry" mergers of elliptical galax- 
ies. In these cases, we collide two remnants of previous disk- 
disk mergers, in order to explore how their properties are 
modified through re-merging. We typically merge two iden- 
tical remnants (i.e. two identical copies of the remnant of a 
given disk-disk merger), but have also explored re-mergers 
of various mass ratios (from 1:1 to 4:1), and mixed morphol- 
ogy re-mergers (i.e. merging an elliptical remnant with an un- 
merged gas-rich disk). In the former case, we generally find 
a similar division in mass ratio at which a "major" merger is 
significant. In the latter, we find the properties are more akin 
to those of other gas-rich (disk-disk) mergers, and the rem- 
nant should for most purposes still be considered the direct 
product of a gas-rich merger. In our re-merger series, we vary 
the orbital parameters, both of the initial gas-rich merger and 
re-merger, and consider systems with a range of initial gas 
fractions in the (pre gas-rich merger) progenitor disks. Our 
re-mergers span a similar range in virial velocities and final 
stellar masses to our gas-rich mergers. 

Each simulation is evolved until the merger is complete 
and the remnants are fully relaxed, then analyzed following 
ICox etal.ld20Q6bh in a manner designed to mirror the meth- 
ods typically used by observers (for details see Paper II). The 
projected surface brightness profile is fitted with stand ard el- 
liptical isoph otal fitting algorithms (e.g. IBender et al.lll987l: 
iBenderl 1988ft and radial deviations of the iso-density contours 
from the fitted ellipses determine the boxyness or diskyness of 
each contour (the 04 parameter). Throughout, we show pro- 
files and quote our results in terms of the major axis radius. 
The effective radius Rf. is the projected half-mass stellar ef- 
fective radius^, and the stellar mass refers to the total stel- 
lar mass of the galaxy. When we plot quantities such as Rg, 
we typically show just the median value for each simulation 
across ~ 100 sightlines. The sightline-to-sightline variation 
in these quantities is typically smaller than the simulation-to- 
simulation scatter, but we explicitly note where it is large. 

3. THE DATA 

We compare our simulations to and test our predictions on 
an ensembl e of observed s urface brightness profiles of ellip- 
ticals from iKormendv et~a l. (2008) and Laueret al. (2007a). 
The important aspects of the observations for our purposes are 
also summarized i n Paper II. Briefly, first is the V-band Virgo 
elliptical survey of lKormendv et al.T(l2008h . based on the com- 
plete sample of Virgo galaxies down to extremely faint sys- 
tems in Binggeli et al. ( 198 5) (the same sample studied in 
ICdte et aLl l2006t iFerrarese et aP I2OO6O . The HST images 
alone, while providing information on the central regions, typ- 
ically extend to only ^ 1 kpc outer radii, which is insufficient 

^ The results here are primarily from equal-mass mergers; however the 
behavior does not change dramatically for mass ratios to about 3:1 or 4:1, 
appropriate mass ratios for comparison with the observations of ellipticals 
used in this paper. At higher mas s ratios, the result is a small bulge in a 
still disk-dominated galaxy (see e.g. lYounger et alJ2008l : IHoDkins et al.l2009l 
'20081), which we do not study here. 

This differs from what is sometimes adopted in the literature, where Sc 
is determined from the best-fitting Sersic profile, but because we are fitting 
Sersic profiles to the observed systems we usually quote both. 



Dissipational Signatures in Core Ellipticals 



5 



to fit the outer profile. iKormendv et alj ( l2008l) therefore com- 
bine observations from a large numb er of source s (including 
Bender et al."1988l "2007"; 'Caon etal."1990', ^199^, 'Dav is et all 
1985; Kormendv et al. 2005; Lauer 1985a; Lauer et al.l l 19951 
20051: iLiu et al. 2005; Pe letier et al . 1990) and new photom- 



etry from McDonald Observatory, the HST archive, and the 
SDSS for each of their objects which enables accurate surface 
brightness measurements over a wide dynamic range: profiles 
spanning ^ 12-15 magnitudes in surface brightness, corre- 
sponding to a range of nearly four orders of magnitude in 
physical radii from ~ lOpc to ^ lOOkpc. 

We a lso add surface brightness prof iles from iLauer et alj 
(I2007ah. further supplemented by Bender et all (119881) . 
iLauer et al] (l2007ah compile V-band measurements of a large 
number of nearby systems for which HST ima ging of the 
galactic nuclei is available. These include the iLauer et aH 
G005) WFPC2 data-set, the lLaine etai](l2003h WFPC2 BCG 
sample (in which the objec ts are spe cifically selected as 
brigh test cluster galax ies f romlPostman & Lauer ( 1995)), and 
the lLauer et all(ll995h and lFaber et al.l([T997.) WFPCl compi- 
lations (note that issues of completeness and environment are 
not important for any of our conclusions). This extends our 
sampling of the high-mass end of the mas s function, but at the 
cost of some dynamic range in the data: iLauer et aTl (l2007ah 
combine these data with ground-based measurements to con- 
struct profiles that typically span physical radii from lOpc 
to ~ 10-20kpc. W e have compared these data with addi - 
tional profiles used in lBender et all (1 1 988i [1991 1 1 993L 1 1 994 . 
and in some cases subsequently updated. These are more lim- 
ited, extending from ^ 30-50pc out to ~ a few kpc, but they 
allow us to construct multicolor (V, R, I) surface brightness, 
ellipticity, and 04/0 profiles, which we use to estimate our 
sensitivity to the observed waveband and photometric qual- 
ity/dynamic range. 

When comparing to cusp populations and gas-rich merg- 
ers from Paper I and Paper II, we occasionally refer to the 
sample of local rernnants of recent gas-rich mergers from 
iRothberg & JosephI ( |2004|) . Details of these observations are 
summarized in Paper I; they are ^T-band profiles of rem- 
nants spanning a range in merger stage, from ULIRGs and (a 
few) unrelaxed systems to shell ellipticals. As demonstrated 
therein, these systems will almost all become (or already are, 
depending on the classification scheme used) typical ^ L» 
cusp ellipticals. 

Here, we are specifically interested in testing the hypothesis 
that core ellipticals retain some memory of extra light, even if 
they have subsequently been modified by dry re-mergers. We 
therefore divide our sample into those systems which are con- 
firmed via HST observations as being either cusp or core el- 
lipticals. We will specifically focus on the core ellipticals, but 
will compare their properties to those of cusp ellipticals and 
gas-rich merger remnants which we study in detail in Paper I 
and Paper II. We exclude spheroidals, as they are n ot believed 
to form in major mergers as are ell ipticals (e.g. iKormendyl 
1 1 985bl \mM IKormendv et al.ll200"8l) (they also predominate 
as satellites at extremely low masses, and not as core galax- 
ies). We also exclude SO galaxies: these likely form a con- 
tinuous family with low-luminosity cusp ellipticals (and their 
properties are consistent with our gas-rich merger remnants 
in Paper II), but robustly separating the dissipational extra 
light and violently relaxed components in such objects would 
require ideal subtraction of the large-scale disk (and fitting 
three-component models involves large degeneracies). Fur- 
thermore, SO galaxies are almost uniformly classified as cusp 



ellipticals in the literature. 

This yields a final sample of « 110 unique core elliptical 
galaxies spanning masses < 0.1 M» to > lOM,, with fitted 
parameters presented in Table [T] (with a comparison sample 
of « 80 cusp ellipticals and « 50 gas-rich merger remnants 
from Paper I and Paper II). There is overlap in the samples 
used; we have ^ 300 surface brightness profiles for our col- 
lection of unique ellipticals, including (for many objects) re- 
peated measurements in multiple bands and with various in- 
struments. This provides a useful means to quantify error esti- 
mates in fits to these profiles and check for systematic sources 
of error, as the variations between fit parameters derived from 
observations in different bands or made using different instru- 
ments are usually much larger than the formal statistical er- 
rors in the fits to a single profile. For sources with multiple 
independent observations, we define error bars (Table [U for 
each fit parameter representing the ~ 1 cr range in parameters 
derived from various observations, typically from three differ- 
ent surface brightness profiles but in some cases from as many 
as w 5-6 sources (where there are just 2 sources, the "error" 
is simply the range between the two fits)''. 

We convert the observations to physical units given our 
adopted cosmology and compile global parameters and cal- 
culate stellar masses as described in Paper II (our conclusions 
are not sensitive to either of these procedures). The data often 
cover a dynamic range and have resolution comparable to our 
simulations, provided we do not heavily weight the very cen- 
tral (< 30pc) regions of HST nuclear profiles. Experimenting 
with different smoothings and imposed dynamic range lim- 
its, we find it is unlikely that resolution or seeing differences 
will substantially bias our comparisons. They can introduce 
larger scatter, however: the robustness of our results increases 
considerably as the dynamic range of the observed profiles is 
increased. 

Throughout, we will usually refer interchangeably to the 
observed surface brightness profiles in the given bands and 
the surface stellar mass density profile. In Paper I and Paper 
II we test this assumption in simulations, and show that once 
the system is relaxed, the optical bands become good proxies 
for the stellar mass distribution, with < 20% variation in our 
M/L over the entire fitted range of radii, in good agreement 
with what is inferred using the observed color gradients or 
full resolved stellar population data to determine M/L versus 
radius (systems have been corrected or excluded for the ef- 
fects of dust, the observed color gradients are weak, age and 
metallicity gradients tend to offset in M/L, and the ages of 
the observed systems are uniformly old, all leading to rela- 
tively little bias). Furthermore, comparison of systems ob- 
served in different bands demonstrates that our conclusions 
are unchanged regardless of the observed bands in which we 
analyze these systems. 

4. LIGHT PROFILE EVOLUTION IN RE-MERGERS 

4.1. Physical Consequences of the Re-Merger 

We first examine how the surface brightness profiles of typ- 
ical merger remnants are modified by re-mergers. To lowest 
order, since the gas is mostly exhausted, the re-merger will be 
dissipationless, and merging two dissipationless, similar sys- 
tems Ml and M2 on a parabolic orbit should roughly preserve 

' In many cases the different observations are comparable; in some there 
are clearly measurements with larger dynamic range and better resolution: the 
errors derived in this manner should in such cases be thought of as the typical 
uncertainties introduced by lower dynamic range or less accurate photometry. 
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their profiles. This leads to the energy conservatio n equation 
(iHausman & OstrikeiifT978l: iHernquist et al.lll993h 

Ef = k{Mi +M2)<j}=Ei = kMi a\ + kM2al (1) 

where cr/ is the velocity dispersion of the final remnant, and 
A: is a constant that depends on the shape of the profile. For 
a 1 : 1 merger of initially similar systems, then, the mass dou- 
bles and the velocity dispersion is unchanged, ct/ = cr,. This 
implies that the effective radius must double as well. Mod- 
ulo a simple overall rescaling of the profile, we therefore have 
reason to believe that our two component approach, applied 
successfully to cusp ellipticals in Paper II, should be able to 
identify the original dissipational and dissipationless compo- 
nents of the pre-re-merger galaxies. 

Figure [T] compares this simple expectation to a subset of 
our re-merger simulations. We plot the "final" specific bind- 
ing energy (binding energy per unit mass efjnai of a given stel- 
lar particle after the system relaxes following a 1 : 1 spheroid- 
spheroid re-merger) as a function of "initial" specific binding 
energy (specific binding energy einitiai before the re-merger, 
but after the initial spheroid-forming merger). Equation ^ 
generalized to the individual stars implies that this should be 
conserved in a re-merger, and indeed, we see a tight corre- 
lation of the form eunai ~ Cinitiai- It has been well-established 
for some time that particles tend to preserve their rank order 
in binding energy in dissipationless mergers (Barnes 198^, 
and we confirm this here (with a scatter in eflnai/einitiai of 
^ 0.15 dex). In detail, fitting the correlations in our simu- 
lations gives gfinai ~ 1.1 einitial- The highly bound particles are 
slightly more bound after the re-merger, because we are con- 
sidering the stellar mass here, not the total mass - relatively 
more of the energy is transferred to the less bound outer halo 
compared to the tightly bound inner spheroid (something like 
this must happen, since some material is scattered to posi- 
tive energies in the re-merger, removing energy from the more 
bound particles). 

If the particle orbital distributions are conserved, then, we 
expect from our estimate above that their radii (or, for non- 
circular orbits, their apocentric radii) should approximately 
double after a re-merger. We therefore compare the "final" 
and "initial" radii of stellar particles in the same manner Note 
that although technically an instantaneous three-dimensional 
stellar particle radius has no inherent physical meaning for 
non-circular orbits, even an instantaneous snapshot (e.g. view- 
ing a galaxy at some moment) will capture most particles near 
their apocenter (in the mean over any ensemble of particles, 
effectively one will recover the apocentric radii, because this 
is where any given orbit spends most of its time). We there- 
fore show in Figure [T] the radii at a given instant before and 
after the re-merger but find that measuring a time-averaged ra- 
dius or the apocentric radii averaged over some time interval 
yields nearly identical results. We consider a narrow range of 
initial radii and examine where the stars at these radii end up 
after a re-merger, plotting the median and ± 1 cr range of their 
final radii (each plotted bin contains at least > 100 particles, 
and most contain > 1000). To lowest order, the median final 
effective radii are roughly twice the initial radii. 

In detail, fitting the points shown yields a mean increase 
in radius of a factor of 1.8 (slightly less than 2, because of 
the fact noted above that energy is preferentially transferred 
to less bound material in the outer halo; in fact this factor 1 .8 
corresponds exactly to the eunai =1.1 einitiai change in specific 
energies, given the virial theorem). At small radii, there ap- 
pears to be a deviation from this, but this can be explained 
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Fig. 1 . — The impact of re-mergers on the stellar light profile. Top: Final 
specific binding energy (egnal) of stellar particles after a spheroid-spheroid 
re-merger is plotted relative to their specific binding energy after the first, 
spheroid forming gas-rich merger (tinitiai ; both e are in units of the character- 
istic energy GMj /Rd, where M, and Rf , are the galaxy mass and effective 
radius before any mergers). Each plotted point is the median value for > 100 
particles within ±0.025 dex in einitial, with error bars showing the ±1 o" range 
of efinai. Different colors and symbols denote different simulations. Dotted 
fine represents exact conservation of specific binding energy; solid line is the 
best-fit constant rescaling ejinai ~ 1 ■ 1 t\mti:A (representing the fact that energy 
is preferentially redistributed to the less bound parts of the halo). Middle: 
Same, but for the radii of each stellar particle before (rinitiai) and after (rfj^a]) 
a re-merger (reflecting their pre and post-merger orbital apocentiic radii). 
Solid line represents the expectation from the eflnal ~ ^initial relation and the 
virial theorem: r^nal ~ 1-8 ''initial- Dashed and dotted lines indicate the ±1 o" 
(0.4 dex) and ±2a range. Bottom: Comparison of initial progenitor pro- 
jected surface mass density profile, and that profile with all radii increased by 
a uniform factor of 2 (the total mass also doubled) with that of the re-merger 
remnant (two representative examples shown). 
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entirely by a combination of resolution effects and the fact 
that the distribution is one-sided (there are only positive radii). 
To lowest order, then, the mass profile is preserved, with all 
radii "inflated" or "puffed up" by a factor ~ 2. We show this 
explicitly in the figure; comparing the final profiles with the 
progenitor profiles expanded by this factor Correspondingly 
(given mass conservation), the effective surface brightness /g 
decreases by a factor ~ 2 and a is conserved; however, typical 
elliptical profiles are sufficiently steep that this "puffing up" 
leads to a final surface brightness profile that is everywhere 
(at each fixed radius) brighter than the progenitors (with the 
effect larger at large radii). 

In short, we have recovered the conventional wisdom that 
that particles tend to preserve the ir rank order i n binding en- 
ergy in dissipationless mergers dBarnesI [19881) . which also 
leads to their mass profiles and central mass concentrations 
being preserved, at least in mergers of systems with certai n 
dist ributions, such as the NFW profile ^Navarro et al."1996'), 
the iHernquistI (119901) profile, or th e Vaucouleurs (1948) 
r'/4 law (e.g. lBovlan-Kolchin e t al.' 2005). 

However, for a given initial orbital apocentric radius or spe- 
cific binding energy, there is in fact a wide range of final or- 
bital apocentric radii, which we will show in § [TO]has impor- 
tant consequences for the details of the observed profile. In 
detail, the distribution of final radii for a given initial radius is 
close to log-normal (out to 3 cr in the wings), with a median 
a factor of ^ 1.8 larger than the initial radius and a 1 ct dis- 
persion of 0.38 dex (factor sa 2.5). Most of this scatter owes 
to mixing in specific binding energy (the scatter in efinai/einitiai 
of ^ 0.15 dex translates to a scatter in rgnai/ initial of 0.3 dex), 
and would be present even if individual orbital structures were 
perfectly conserved; the rest owes to mixing altering the na- 
ture of individual orbits (i.e. changing the isotropy of the sys- 
tem and moving stars between tube and box orbits); only a 
small component of the scatter is attributable to the inherent 
uncertainty in defining an orbital "radius" in a non-trivial po- 
tential. Figure [1] shows the 1 a range for each bin in initial 
radius and this fit to the entire distribution. 

Figures |2] & [3] show how the surface density profile of the 
final galaxy is modified by a re-merger. We show the sur- 
face brightness profiles for two gas-rich merger remnants, 
one in which the "extra light" signature in the center of the 
galaxy is plainly obvious (a good match to e.g. NGC 4459 in 
Virgo), and one in which it is more subtle (i.e. the extra light 
blends more smoothly into the outer profile; analogous to e.g. 
NGC 4886a). For each, we also show the profile of those 
stars specifically formed in the central starburst and those not. 
We then show the surface density profiles after a re-merger 
(for each, we show the profile after two different re-mergers 
with different orbital parameters - the choice of orbital pa- 
rameters has only a weak impact on the profile). We can fol- 
low the stars which originally formed in the gas-rich merger- 
induced starburst, and plot their profile as well (even though 
there is no new starburst in the re-merger, since the systems 
are gas-depleted). In both cases, the profiles are roughly pre- 
served, as expected based on the arguments above. The star- 
burst light profile, while stretched by a similar factor 2 by 
the re-merger, remains relatively centrally concentrated and 
continues to dominate the profile at radii < 1 kpc, while the 
non-starburst stars that were originally relaxed into an outer 
Sersic-like profile continue to follow a similar distribution. 

The light profile is, however, smoothed out by this re- 
merger. In the case where the original remnant shows a sharp 



departure from the outer profile, the departure remains rela- 
tively prominent and detectable even after a re-merger But, 
when the original starburst light is more smoothly distributed, 
the re-merged light profile can be smooth, without any ob- 
vious feature. This case demonstrates that because the in- 
ner parts of the outer/non-starburst component shape (already 
being "hot") are not affected much in the re-merger, while 
the starburst component is heated, we often find the starburst 
component is less dominant in the center after the re-merger 
(they have mixed). It is still important to the central profile 
and much more concentrated than the outer light, but the dis- 
tinction is less sharp. It is therefore interesting to test whether 
our empirical fitting routines reliably recover the mass pro- 
file of what was, originally, the starburst light, even after a 
re-merger 

4.2. Decomposing the Remnant Profile 

We therefore adopt a two-component approach, described 
in Paper I and Paper II, to fit the final (post re-merger) surface 
brightness profile of each system. The methodology arrived 
at and tested in these papers amounts to a simple prescription: 
we fit the total observed light profile (temporarily ignoring the 
known true physical decomposition of the light profile) to the 
sum of two Sersic components. Based on the arguments and 
tests in Paper I and Paper II, we fit to an outer Sersic plus 
cusp or extra light model, with an outer component for the 
original pre-starburst/dissipationless stars with a free Sersic 
index, and an inner component reflecting the remnant of the 
dissipational/starburst population. Where data are limited at 
small radii, or in the case of our simulations where we do not 
resolve the structure on very small nuclear scales, we find the 
best results (in an average sense) adopting a fixed Sersic index 

= 1 for the inner component, but we note (as demonstrated 
in Paper II), that the same mean results are recovered for a free 
inner Sersic index, and that where the data are of sufficient 
quality (i.e. for systems with HST nuclear profiles extending 
to very small radii), it is possible to free this parameter. 

Our fitting procedure then takes the following form: 

/tot=/' exp|-/7;,(-^)'''"'| 

^ "extra ^ > 

+/oexp|-/7„(-^)'^"'|, (2) 



where 7?extia and 7?outer are the effective radii of the inner 
(«j ^1) and outer (free Wj) components (which we identify 
with the starburst and old bulge or pre-starburst components, 
respectively), /' and 4 are the corresponding normalizations, 
is the Sersic index of the inner (extra light) component 
(fixed «s = 1 where resolution limits apply) and is the Sersic 
index of the outer bulge or pre-starburst component. The con- 
stant h„ is the appropriate function of such that /?extia and 
^outer correspond to to the projected half-mass radii. We stress 
that these fits are intended to recover the appropriate profile 
and fit where these components are physically relevant - i.e. 
at radii > 30-50pc where our simulations are well-resolved 
(note that the typical effective radii of even the inner compo- 
nent are much larger than this, ^0.5-1 kpc). For most of the 
galaxy observations considered here, this corresponds to radii 
> 0.5", a factor ^ 10 larger than the HST diffraction limit. 

The results of applying these fits are shown in Figure|2] De- 
spite the smearing of the transition to the extra light that owes 
to the re-merger, this approach still reliably recovers both the 
true outer Sersic index of the violently relaxed non-starburst 
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Fig. 2. — Top: Surface mass density of a gas-rich merger remnant (black), physically decomposed into stars formed prior to the final merger (which are then 
violently relaxed; red) and stars formed in the dissipational starburst (blue). The pre-starburst component follows a Sersic law closely, with best-fit index shown, 
and the starburst mass fraction is labeled. The second derivative of the light profile (d^/j/(dlogr)^) is shown below. We compare this with our two-component 
(Sersic plus cusp or extra light) fit (inner exponential and outer Sersic) to the total light profile, with the Sersic index of the outer component and mass fraction 
of the inner component, and rms scatter (cr^) about the fit. Residuals from the fit are shown below. Middle: The same, after a spheroid-spheroid re-merger of 
the remnant above. The (post re-merger) profile of the stars originally formed in the starburst is now shown, and the re-merger remnant is fitted with the same 
procedure. Bottom: Same, but for a re-merger with different orbital par'ameters. We show a system with a prominent transition to the extra light (left panels), 
and a more typical system with a smoother transition between outer and extra light (right panels). Our two-component, cusp plus Sersic function fit accurately 
recovers the profile of the violently relaxed component and mass fraction of the starburst component, even after a re-merger The remnant of the original starburst 
(extra light) is less obvious, however, after re-mergers. 

does not intierently rely on a "break" in tlie profile - if there 
is a smooth but significant change in e.g. the effective con- 
centration of the system going from inner to outer regions (or 
equivalently the rate of change of the logarithmic slope), it 
will be reflected in the fits by a dense inner component and 
more extended outer component. 

Figure m shows the result of these fits to our gas-rich disk- 
disk merger and gas-poor re-merger simulations. We directly 
compare the fitted extra light mass fraction to the mass frac- 
tion and size of the known physical starburst component, and 
find that the fitted components recover the appropriate values 
in the mean with a factor ^ 2 scatter. This result is robust with 
respect to e.g. the mass, orbital parameters, mass ratios, ini- 
tial gas content, treatment of feedback and model for the ISM 
equation of state, and redshift of our simulations. We refer 
to Paper I for a more detailed comparison where we similarly 
demonstrate the ability to recover the starburst effective radius 
and outer component Sersic index. There is no significant bias 
evident in our re-merger sample: although it is considerably 
more limited, the representative systems we have re-merged 
appear to lie in a similar location in this space to their progen- 
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Fig. 3. — As Figure |2]( bottom panels); showing the re-merger remnant 
profiles and physical breakdown, in p — r^^'^ instead of p — logr. The smooth 
case (right) retains its extra light in a physical sense (and it is successfully 
recovered by our decomposition procedure), but the profile shows no obvious 
upward "break" relative to an outer Sersic index. 

Stellar material and the mass fraction in the dissipational cen- 
tral starburst. We have tested this for our entire library of 
re-merger simulations, and find that while the smoothing of 
the profile makes the object-to-object errors somewhat larger, 
it does not alter the ability of this method to recover, on aver- 
age, the appropriate breakdown between physical components 
of the galaxy light profile. This is in part because the fitting 
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Fig. 4. — Success of our proposed two-component empirical decomposition at recovering ttie Icnown pfiysical parameters of the galaxy starburst and pre- 
stai'burst (violently relaxed) components. Left: Mass fraction in the fitted "exfi'a light" component /extra versus the known mass fraction of the physical starburst 
/sb, in our gas-rich (disk-disk) merger simulations. Each point is the average across ~ 100 sightlines to a given simulation; median sightline-to-sightline variance 
is shown (fa 0. 15dex; plotted error bar). Various colors and symbols denote different initial disk gas fractions, orbital parameters, and merger redshifts (key). 
Center: Same, but comparing our gas-rich simulations (grey, for clarity) with our gas-poor (dry) re-merger simulations (magenta). The starburst fraction /sb 
represents the sum of the starburst masses in the two ellipticals re-merged. Symbol shapes denote original orbits of the progenitors in the gas-rich mergers 
(although the re-mergers follow a mix of orbits). Each of these points is close it its progenitor gas-rich merger (progenitor initial gas fractions can be inferred 
from the location relative to points in left panel; this is less significant than /sb). Right: Distribution of /cxtra/,/sb> including sightline-to-sightline and simulation- 
to-simulation variation. Histograms are for our gas-rich (grey) and re-merger (magenta) simulations. Dotted blue line shows a Gaussian with cr = 0.27 dex. The 
fitted /ex(ra recovers the physical /sb on average, with a factor ^ 2 scatter, and without any significant bias from any varied simulation parameters. Re-mergers 
yield a similar result - the error bars only significantly increase for systems with either very little starburst components (where the smoothing is significant; see 
Figure |2) or a large number of re-mergers. 



itors, with similar scatter At the very lowest starburst frac- 
tions, the scatter does increase in re-mergers, owing to the 
smoothing of the profile and mixing of the (small) extra light 
mass with the outer dissipationless mass (shown in Figure|2]i. 

We note that these simulations do not directly model the 
process of "scouring," or the gravitational scattering of in- 
dividual stars by a coalescing black hole binary. Unfortu- 
nately, it is prohibitively expensive to model these processes 
self-consistently in a galaxy-scale simulation, as they require 
resolving the masses of individual stars and sub-parsec spa- 
tial scales. We instead adopt a standard gravitational soft- 
ening prescription, and for this reason the predicted profiles 
should not be taken literally within the simulation smooth- 
ing lengths (~ 50pc). In fact, scouring is expected to alter 
the galaxy profile significantly only near the radius of in- 
fluence of the central black hole, i.e. <^ 30-50pc for most 
of the systems of interest (where the enclosed mass of the 
galaxy is ~ Mbh; as generally expected from theoretical mod- 
els and simpl e scaling arguments; Milosavlievic et al. 2002; 
lMerrittll2006l; ISesanaet al .1120071) . which is below our simu- 
lation resolution limits and well below the relevant size scale 
at which the extra or starburst light we model is important 
(^ 0.5 - 1 kpc). While a very small fraction of stars at these 
radii may pass near the black hole on centrophilic orbits, there 
is no scenario, in any plausible scouring calculation, in which 
the profile at radii ~ 0.5- Ikpc, enclosing a stellar mass 
> IOOMbh, could be sufficiently altered by this particular pro- 
cess in order to change our conclusions that the original extra 
light is recovered after re-mergers. In other words, scouring 
flattens ("cores out") the central peak in the dissipational or 



"extra light" profile - it changes the exact slope/shape of the 
central light profile as r ^ 0, but not its total mass or ^kpc 
extent (and we have tested in numerical experiments that our 
adopted fitting procedures are robust to changes in the nuclear 
profile shape). In short, the effects of "scouring" on this spe- 
cific calculation are no different than the effects of our finite 
spatial and mass resolution. 

Figure |5] summarizes the physical properties of the surface 
brightness profiles of gas-rich merger remnants and gas-poor 
("dry") re-merger remnants. The original gas-rich merger 
remnant is fundamentally two-component, with a Sersic-like 
violently relaxed envelope (with characteristic Sersic index 
«j ^ 2.7 ± 0.7, discussed in Paper II) and an inner dissi- 
pational/extra light component that dominates the light pro- 
file inside of ^ 0.5 - 1 kpc. The inwards extrapolation of 
this extra light forms the central cusp, a (roughly) power- 
law like continuation (with characteristic power-law slopes 
dln//dlnr -0.7 ± 0.3). The extra light, and by exten- 
sion the cusp (insofar as it is simply the inwards continu- 
ation of the extra light profile), as well as embedded kine- 
matic subsystems (see ICox et all l2006bl; iNaab et akl l2006al) 
which make the remnant isophotal shapes more disky and 
yield rapid rotation, are formed by the dissipational star for- 
mation event(s). Rapid rotation and disky structure are also 
frequently observed in what we now recognize as the extra 
light regions of cuspy ellipti cals (e.g.. lEmsellem et al.ll2004l 
l2007tlKormendv et al.ll2008h . 

In a dissipationless re-merger, both the original dissipa- 
tional and dissipationless components are "puffed up" by a 
uniform factor ^ 2. The scattering of stars (the scatter in Fig- 
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Fig. 5. — Summary of the physical properties of spheroid profile shapes (here for a typical ~ L, elliptical), as formed in a major, gas-rich merger (left) and 
then modified by a major gas-poor, spheroid-spheroid re-merger (right). Lower panels show the logaiithmic derivative of the total (observed) profile. In the 
original merger, stars from the pre-merger disks are violently relaxed into an extended Sersic-like envelope (n, ~ 2.7). Conservation of phase-space density 
prevents this component alone from reaching the high densities of observed ellipticals. Gas dissipation, however, yields a nuclear starburst, leaving a dense "extra 
light" component on top of the outer component, dominating the profile at r < 0.5 — 1 kpc. The nuclear "cusp" is the inwards continuation of this dissipational 
component. Together this yields a global profile with typical «j ~ 4 and densities of observed ellipticals (much higher than their progenitor spirals). In a re- 
merger, both components are "puffed up" by a factor ~ 2, and scattering of stars broadens the original dissipationless/envelope component (to /jj ~ 3.6). There 
is no new dissipation, but the remnant of the original dissipational component continues to dominate the profile within ~ 1 kpc and is (in principle) separable. 
Black hole "scouring" scatters stars from the nuclear cusp/extra hght - evacuating ~ Mbh worth of stai's and flattening the central profile to form a nuclear "core." 
Although the apparent effects on the extra light profile can be large, the total mass in the "scoured" region here (the "missing mass") is only ~ 2% of the exlrci 
light mass (most of which is near ~ 0.5 - 1 kpc). An'ows denote the typical resolution limits of our simulations (sufficient to model the structure of the extra light 
and outer components, but insufficient to simultaneously resolve the transformation of cusps into cores). 



ure[T]l makes the envelope of the original dissipationless com- 
ponent broader, raising its Sersic index (discussed in §|9]l. The 
mass of the original dissipational component is conserved, as 
is its dominance within the central ^kpc - this is a conse- 
quence of the overall preservation of profile shape. However, 
heating of the inner component and mixing of stars according 
to the scatter in Figure [1] makes the relative prominence (the 
rise above the outer profile) somewhat more smooth. With no 
scouring mechanism, the central slope would still eventually 
rise about as steeply as the progenitor cusp (again, reflecting 
the conservation of profile shape). Scouring, or scattering of 
stars by a nuclear binary black hole, will however evacuate 
the central regions and flatten this slope. 

Here we have applied a toy scouring model following 
lGualandrir & Merritt (2007), which kicks out a total of ^ 
2Mbh worth of stars. The scoured region (i.e. the scoured or 
missing mass) is therefore represented by the area between the 
scoured and un-scoured profiles. Because the extra light dom- 
inates the stellar density in the central regions (and formed 
the original cusp), scouring acts primarily on the stars in the 
nuclear region of the extra light. This leaves a central core 
with a shallow slope (dln//dlnr^-0.2±0.2). The re-merger 
also tends to des troy the embedded kinematic subsystems (see 
ICox et al.ll2009l) and no new source of dissipation is avail- 
able to re-form them, leaving a remnant with boxy isophotal 
shapes and less rotation. 



4.3. Caveats and Limitations 



We emphasize here the clear hierarchy of scales: the large 
scale profile and envelope is made primarily from the dissipa- 
tionless violently relaxed component. The extra light domi- 
nates the profile within '^kpc scales, a factor of a few to an 
order of magnitude smaller than the effective radius. Pro- 
cesses such as scouring and the observational distinctions be- 
tween cusps and cores occur at scales yet another order of 
magnitude (or more) smaller, < lOOpc. Because of this hier- 
archy, although scouring can appear to alter the shape of the 
"extra light" profile significantly, it has a negligible effect on 
the estimate of the total extra light mass fraction (here, the 
entire mass content inside the "scoured" radius is just ~ 2% 
of the extra light mass) or effective radius (this 7?^ ^ 1 kpc). 
In general, we can continue to treat galaxy profiles as multi- 
component (in a physical sense) after a moderate number of 
dry mergers, because the mixing of stars over factors ^ 2 in 
radii in major re-mergers is insufficient to eliminate the dif- 
ferences between such different scales. 

This will no longer be true after a sufficiently large number 
of re-mergers (if there is an initial order-of-magnitude differ- 
ence in the scales of the extra and outer light, then simple 
scalings suggest ^3 1:1 mass ratio re-mergers, or 5- 10 
more likely 1:2 - 1:3 mass ratio re-mergers, could com- 
pletely "blend" the dissipational and dissipationless compo- 
nents). However, cosmological simulations, clustering mea- 
surements, and halo occupation models (Mailer et al. 20061; 
IZheng et ah.2007 ; Mas j edi et al. 2008; Ho pkins et al. 2008d) 
suggest this will only be important in the most extreme M* > 
KV^Mq BCG populations (which constitute only ~ 10-20% 
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of the mass density within the core elliptical population, and 
only ^ 3-5% of the mass density in ellipticals). In this 
regime, there are other reasons to treat our models with cau- 
tion: such systems almost all experienced their first mergers 
at very high redshifts where the progenitors are less well- 
understood, they live in unusual environments, and at these 
masses it may be the case that growth by a large series of mi- 
nor mergers becomes more important than growth by major 
mergers. We do not intend, therefore, for our models here to 
be considered as analogs for these systems; our intention is to 
understand the bulk of the core elliptical population at masses 
< a few L„ (constituting 80-90% of the mass density in 
core elUpticals), the vast m ajority of which stil l live in less 
extrerne environments (e.g . iBlanton et al.1 120051: IWang et al] 
120061: iMasiedi et al J 120061) and are expected to have experi- 
enced a small number (~ a couple) of major re-mergers since 
redshifts z ^ 2-4. 

Ellipticals with cores, then, if they are the products of the 
re-mergers of cuspy progenitors, should have just as much 
"extra light" (dissipational content) as cusp ellipticals, and we 
will show in §|6]and §|2]that this surviving dissipational com- 
ponent obeys the same correlations and has similar proper- 
ties to the extra light in observed cusp systems. However, the 
term "extra light" has not generally been applied to core ellip- 
ticals, because observations typically show less (compared to 
what is seen in cusp ellipticals) of an obvious upward break 
in the profile at small radii relative to the inwards extrapola- 
tion of their outer Sersic profiles. In fact, core ellipticals are 
sometimes referred to as "missing light" ellipticals, because 
they show a deficit in their central light profiles, relative to 
the extrapolation of an outer Sersic profile. These differences 
arise naturally from the smoothing and mixing processes in 
re-mergers, and in fact are completely consistent with our ad- 
vocated formation scenario (see 

As discussed above, "extra light" should be physically 
identified with stellar populations originally formed in dis- 
sipational starbursts (on top of more extended violently 
relaxed populations from dissipationlessly merged stellar 
disks), whereas "missing light" is associated with the deficit 
of stars in the nucleus owing to scattering of some small mass 
by a merging binary black hole. In a physical sense then, el- 
lipticals can be both "extra " and "missing " light ellipticals 
- the terms as we mean them are not mutually exclusive. In- 
deed, if "extra light" means stars from dissipation (on top of 
dissipationlessly assembled stars) - as we use the term - then 
all ellipticals must have "extra light" at some level. And if a 
black hole binary does not merge sufficiently rapidly (if there 
is little gas in a recent major merger), scouring will happen. In 
this physical scenario, cored or "missing light" ellipticals are 
re-merged extra light ellipticals where the inwards extrapola- 
tion of the original extra light (which constituted the original 
"cusp") has been flattened by scouring, forming a core in the 
center of the dissipational extra light profile. 

For these reasons, we will continue to refer to our fitted 
central components as "extra light" and interpret them as the 
surviving dissipational components from the gas-rich mergers 
that formed the core elliptical progenitors. Given our interpre- 
tive context and calibration from our re-merger simulations, 
we will show that such components appear ubiquitous in core 
ellipticals and obey well-defined scaling laws closely related 
to those of the dissipational components in cusp ellipticals. 

However, that is not to say that core ellipticals would nec- 
essarily be identified as "extra light" galaxies in the sense 
of certain traditional observational metrics. Specifically, the 
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Fig. 6. — Illustration of the typical effects of a re-merger on the stellar 
population gradients of spheroids. For four distinct original gas-rich merg- 
ers (a < L, and > L^, case with typical starburst mass fractions /jj, 10%, 
and a highly gas-rich and very gas-poor ~ L, case, as labeled), we show the 
metallicity as a function of radius (plotted line is the median, variation across 
~ 100 sightlines is Ri 0.05-0.1 in [Z/H] at each r). We show the B-band 
light-weighted metallicity. at each of three approximate times after the major 
starburst corresponding to the original gas-rich merger (the gas-rich merger 
remnant is evolved passively, the re-merger case experiences a re-merger be- 
tween 1 = — 3 Gyr). At the same mean stellar population age, the re-merger 
has had little effect on the profile, washing it out only slightly. The starburst 
mass fraction and overall stellar age are much more important to the gradient 
strength. We find similar results for age, [a/Fe], and color gradients. 

smoothing of the transition between inner and outer com- 
ponents in re-mergers and raising the outer Sersic profile n^ 
by scattering stars to large radii means that the extra light 
may appear less pronounced. Re-mergers also tend to de- 
stroy kinematic subsystems and mix stellar orbits (Cox et al] 
|2009|) , potentially wiping out obvious changes in isophotal 
shapes or kinematics near the radius where the estimated ex- 
tra light begins to dominate the profile. Absent such features 
or a physical model (as we have in the form of our simula- 
tions) to lend an interpretive physical context and suggest a 
multi-component nature in the first place, there will always 
be an infinite space of functional forms that fit the observed 
profiles with arbitrarily high accuracy and do not include 
an explicit "extra" or secondary component. This does not 
mean, of course, that the remnant dissipational populations 
are not manifest in the fits (after all, both our fit and others 
equally represent the data to a meaningful physical accuracy 
"■^i ^ 0.1 magarcsec"^, comparable to the inherent point-to- 
point variance in simulations and observed systems). The 
parameters of interest here will simply be reflected in more 
indirect fashion (in some combination(s)) in whatever fitting 
function is adopted. We discuss some of these comparisons 
with alternative fitting functions commonly applied to core 
ellipticals in the literature in §|9] Given our aims and desire to 
compare simulations and observations on a uniform footing, 
however, we proceed with this specific machinery as a means 
to interpret the observed profiles of core ellipticals. 

5 . IMPACT ON STELLAR POPULATION GRADIENTS 

In Paper II, we demonstrated that dissipation can give rise 
to strong stellar population gradients in ellipticals, and that 
this relates e.g. the gradient strengths, sizes, and extra light 
masses. We therefore briefly examine how these gradients 
might be modified in re-mergers; and show that the gradients 
tend to survive, preserving our predictions from Paper II. 

Figure |6] shows the metallicity as a function of radius from 
several representative simulations, comparing the evolution 
with time of gradients in simulated ellipticals that evolve pas- 



12 



Hopkins et al. 



_ 0.3 



* - 


* ♦ ^ * 

1* * ° « * f J 


* 

♦ 

° ^^ o : 


• Gas-Rich Mergers 
r *^ <> '-^ Re-Mergers 

L * o * 



0.0 0.1 0.2 0.3 0.4 1 10 

[ZJH] Age„ [Gyr] 



Fig. 7. — Stellar population gradients in metallicity or age as a function 
of central metallicity or age, in a subset of gas-rich merger remnants (filled 
symbols) and re-mergers of these same systems (open symbols). Color de- 
notes the initial gas fraction of the original gas-rich mergers, and symbol 
type denotes orbital parameters, as in Figure |4] The dilution of gradients in 
re-mergers is weak - they continue to trace the same trends as gas-rich merg- 
ers of the same stellar age. The effect of increased dissipation in the original 
merger is dominant, and the re-merger is less important in general than the 
simple subsequent aging of the stellar populations by ~ a few Gyr after the 
gas-rich merger. 

sively after the original gas-rich merger to (originally identi- 
cal) systems that have had another (dry) re-merger in the in- 
tervening time. The details of our methodology are described 
in Paper 11, but, briefly, we model every stellar particle as a 
single burst population with an age and metallicity given self- 
consistently by the star formation model in the simulation. 
Stars formed before the original gas-rich merger (i.e. those in 
our initial conditions) are placed on the mass-metallicity re- 
lation for disks of the appropriate mass at the corresponding 
redshift for our simulation. We project the remnant (along 
~ 100 sightlines, but the results depend only weakly on view- 
ing angle) and calculate the light-weighted metallicity, stellar 
population age, and other properties in small radial annuli. 

Naively comparing the profile in e.g. metallicity or age after 
a re-merger to that immediately after a gas-rich merger, it ap- 
pears that there is a significant difference. However, this owes 
largely to age effects (the stars in the re-merger have aged sig- 
nificantly during the merger). Clearly the correct procedure is 
to compare at fixed age or post-original gas-rich merger time. 
It is also important to rescale the systems by their effective 
radii, as we expect a roughly uniform puffing up of the rem- 
nant in a re-merger 

Figure |6] shows that when we do this, there is relatively lit- 
tle smearing out of the gradients by re-mergers. There is some 
effect, but it is generally of the order of the scatter in the gra- 
dients across sightlines, and is generally weaker than the ef- 
fects of age or variation in initial gradient strength across our 
simulations owing to different degrees of dissipation. This is 
expected - Figure [T] demonstrates that the rank order in radii 
is preserved in a mean sense. Although there is some scatter 
in the particle radii, it is not dramatic as far as gradients are 
concerned. Note in Figure |6] the change in metallicity with 
radius is quite smooth, typically decreasing from the central 
metallicity over ~ 2 orders of magnitude in radius. Compared 
to this dynamic range, the typical dispersion in final radii for 
stellar particles at a given initial radius (~ 0.4dex) is signifi- 
cant, but not large. Thus even strong gradients are relatively 
robust to the smoothing effect in re-mergers. 

Figure|7]summarizes these results for our re-merger simula- 
tions, plotting the gradient strength (quantified as the logarith- 



mic slope in metallicity or age versus radius at ^ 0.1 - IR^; 
for details see Paper 11) as a function of central (averaged 
within an aperture of /?f./8, comparable to typical observa- 
tions) metallicity or stellar population age. We show both 
gas-rich merger remnants and the products of dry re-mergers 
of those same systems. Quantified in this manner, there is 
no significant offset between the re-mergers and the passively 
evolved gas-rich merger remnants. 

In Paper 11 we showed that these distributions and the 
gradient strength as a function of mass and velocity disper- 
sion, as well as color gradients as a function of elliptical 
age, all agree well with those observed in elliptical popu- 
lations. Because these are not entirely washed out in re- 
mergers, this agreement continues to hold for re-mergers and 
core ellipticals. It has been specifically noted observati onally 
(e.g. Lauer et al. 2 005; Ferrare se et al. 2006; Kuntsch ner et alj 
2006; McDermidetalJ |2006|) that core ellipticals and even 
BCGs exhibit typical color and stellar population gradients 
compared to other galaxies of similar mass and age. To the 
extent that they may represent re-mergers of (originally) cusp 
ellipticals, our simulations imply that this is the natural ex- 
pectation. 

6. DISSIPATIONAL CONTENT IN ELLIPTICAL GALAXIES WITH 
CORES 

6.1. Fitting Surface Brightness Profiles 

If the preceding discussion is correct, extra or starburst light 
survives a galaxy merger in some sense. We expect, then, 
that even if core ellipticals have been shaped by spheroid- 
spheroid re-mergers, and their profiles smoothed, they should 
show a central stellar density similar to our simulations cor- 
responding to what was (originally) their starburst light con- 
tent, which should be recovered (on average) by our applied 
decompositions. Are observed systems consistent with this 
scenario? 

Figure[8]shows the surface brightness profiles of three core 
ellipticals from moderate mass (NGC 3379, a ~ core el- 
liptical) to very massive (NGC 4874, a > lOM, elliptical). In 
these systems, there is a distinct shoulder in the light profile, 
which is evident in a completely non-parametric sense in the 
first and second logarithmic derivatives of the profiles - the 
slope of the profile is clearly non-monotonic with radius over 
a brief range. Moreover, the shape of this feature is similar 
in all three cases. There are perhaps alternative explanations 
for these features (NGC 3379 may in fact be a face-on SO, see 
De Lorenzi 2008, and NGC 4874 is a very massive galaxy 
with a large envelope, albeit not technically a BCG) but these 
are not unique systems - there are a significant number of 
other such cases, as we describe below. 

The features in these observed light profiles are similar, 
both in their shape and characteristic radii, to the physical 
transition in our simulations where the extra or starburst light 
begins to dominate with respect to the dissipationless, non- 
starburst light. Applying our two-component fit decomposi- 
tion to these profiles, we find a good fit that (as expected) is 
able to explain these shoulder-like features as the transition 
between extra/starburst light and the outer light distribution. 
In a parametric sense, if we fit a Sersic or Nuker-law profile 
to the central regions (specifically, we consider radii between 
lOOpc - within which the core structure is important - and 
~ 0.5 - 2 times the radius at which our two component fit de- 
termines that the two components are equal; the exact range 
is not important), we find that the substantial curvature of the 
inner profile leads to a fit which falls short of the shoulder. 
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Fig. 8. — Examples of core ellipticals with possible surviving signatures of starburst or "extra" light from cusp progenitors in their nuclei (preserved, albeit 
smoothed out, from the original, gas-rich spheroid forming merger). Top: Observed surface brightness profile (open circles) and best-fit two-component model 
(solid line). The inner (starburst remnant; black dashed line) and outer (violently relaxed dissipationless remnant; black dotted line) components of the best-fit 
model are shown, with the radius (vertical dotted hne) at which the two are equal. Nuker fits (red dot-dashed) to the profile between the central ~ lOOpc (inside 
which the core dominates the fit) and radius of the shoulder are shown - these fall short where our fitted outer component dominates the light, reflecting a 
significant change in the curvature of the profile. Middle: First derivative of the surface brightness with respect to logr. Bottom: Second derivative of the profile. 
Characteristic features or shoulders in the profiles can survive re-mergers, manifest plainly in a non-parametric fashion in both first and second derivatives, 
conesponding to the radius at which our method estimates the starburst component begins to dominate. (The specific objects are discussed further in the text.) 
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Fig. 9. — As Figure[8] but for systems with less prominent/obvious features. Despite being weaker, a similar transition is seen, and the amplitude in the first 
and second profile derivatives is comparable. The Nuker fits to the central profile similarly underpredict the outer component, as expected if a second physical 
component becomes important. 
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and fails to reproduce the outer profile dominated (in our two- 
component model) by an outer, violently relaxed population. 

Figure |9] compares three galaxies with relatively weak fea- 
tures, in Virgo. These are systems where presumably the scat- 
tering of stars in re-mergers, discussed in § |4] has smoothed 
any strong break to the extra light that might have been promi- 
nent in the progenitors. Despite not being obvious on vi- 
sual inspection, however, the first and second derivatives of 
the profile show similar signatures, and in fact the deviations 
from a single profile fit (a Sersic or Nuker law) extrapolated 
outward from the central region (interior to the shoulder) are 
significant. Moreover, we find here (and in Figures |2]|5]) that 
although the formal uncertainties may be larger when the pro- 
files are smoothed in this manner, our two-component model 
is able to successfully fit them: § |4] demonstrates that even 
in systems with smoother profiles than these, our applied de- 
compositions are robust. 

Together with the arguments in § |4l this gives us reason- 
able confidence that we can extend our analysis to the com- 
pilation of observed ellipticals described in § [3] Figures [TOl - 
[T4l show surface bright ness profiles of objects in the Virgo 
core elliptical sample of iKormendv et al] (l2008l) . in order of 
most massive to least massive^". For each object, we plot 
the surface brightness profile with the best-fit two compo- 
nent model, and the corresponding fitted outer Sersic index 
and extra light fraction. We refer to Paper 1 for the same 
comp arisons with local observ ed gas-rich merger remnants 
from iRothberg & JosephI (|2004|) . and Paper 1 1 for the compar- 
isons w ith the cusp ellipticals in Virgo from iKormendy et al.l 
((2008"). A complete list of fit parameters and compiled galaxy 
properties is included in Table [T] 

We also compare each observed system with our library of 
re-merger simulations, in a non-parametric fashion. We do 
this by allowing the normalization of the simulated light pro- 
file to vary (within ±0.5 dex), and quantifying the (vari- 
ance of the observed points with respect to the simulated light 
curve at > 1 gravitational softening length) of each simula- 
tion. We allow the normalization to vary because we have a 
finite number of simulations and therefore do not sample a 
continuum in e.g. total brightness, but instead discretely sam- 
ple at factor ^ 2 intervals (we do not allow the simulated pro- 
files to vary by more than this amount, to avoid an unphysical 
match to a simulation with very different total mass). We do 
not allow any other parameters to vary - i.e. we allow limited 
rescaling in the surface brightness of the simulated galaxies, 
but not their radii or other properties. Despite the allowed 
surface brightness rescaling, the best-fit simulations almost 
always have similar total luminosities to the observed system, 
because they must have a similar effective radius in order to 
be a good match. Considering 100 sightlines to each of 
our simulations (although, as noted in Paper 1, the observed 
surface brightness profile varies by a small amount sightline- 
to-sightline), we find the best fit to each observed system. 

We show in Figures [T0lfT4l the three simulations which most 
closely match the observed light profile. For the best-fit sim- 
ulation, we also show the profile of the stars formed in the 
original, central, merger-driven starburst, as described in § [T] 
We show in the figures the outer Sersic indices fitted to these 
simulations, along with the typical range both across sight- 
lines and across the best-fitting simulations (which together 

NGC 4261 is technically in the Virgo W cloud , but we s how i t be- 
cause it is included in the Virgo samples from .Ferrarese et alj j2006l) and 
IKormendv etaH <2008h . 



give some rough approximation to the range of n , which might 
be observed for these galaxies along different sightlines). We 
also show the best-fit starburst mass fraction, along with the 
range across the best-fitting simulations (described below), 
and the rms residuals of the observed points with respect to 
the best fit''. In nearly every case, we easily find simula- 
tions which provide an excellent match to the observed pro- 
files, with variance often less than even a multi-component 
parameterized fit. The fits are good over the entire dynamic 
range from the largest observed radii 100 kpc) down to our 
resolution limits 50 pc)'^. 

In addition, for these simulations we show the isophotal 
shape and ellipticity as a function of major axis radius, com- 
pared to that observed. Note that we do not fit these quantities, 
only the surface brightness profile. We show, for each simu- 
lation, the range across sightlines in these quantities - it is 
clear that these depend much more strongly on sightline than 
the surface brightness profile (this is primarily why we do not 
fit these quantities). In every case, there is a significant frac- 
tion of sightlines with shape and ellipticity profiles roughly 
consistent with those observed, but the simulations highlight 
the range of profile shapes for similar spheroids to those ob- 
served. 

Our simulations are also consistent with the observed kine- 
matic (V, a) profiles of these systems, but we do not show an 
explicit comparison for individual objects because sightline- 
to-sightline variations are sufficiently large (see § [D as to 
yield only weak constraints. For more robust quantities, such 
as e.g . the weighted rotation Ar defined in lEmsellem et alJ 
(l2007h . we find that our re-merger remnants yield good agree- 
ment with the (nearly universal) radial Xr{R) profile of "slow 
rotators" in the SAURON sample (see their Figure 2). The 
distribution of global rotation properties (e.g. V /a) and corre- 
lations with isophotal shapes can be constraining (for a popu- 
lation): we consider this in §[8] There is, however, a great deal 
more kinematic information than just a radial mean velocity 
profile: more detailed constraints for individual objects can, in 
principle, be obtained by comparison of simulations and high- 
resolution two dimensional velocity field maps (with higher- 
order velocity moments such as h^, and /i4). Such a compari- 
son is clearly warranted, but is outsi de the scope of t his paper 
and will be the topic of future work jCox et alj|2009l) . 

Figures [TSlfTT] again show the observed and best-fit sim- 
ulated surface brightness profil es, for a subset of c onfirmed 
core ellipticals in the sample of iLauer et al.l (l2007a !). in order 
of V-band magnitude from brightest to faintest. For each sim- 
ulation in our library, we have a corresponding to its good- 
ness of fit to the observed profile, and the genuine physical 
starburst mass fraction /sb from the original gas-rich merger. 
We can therefore construct a x^-weighted distribution of /sb 
for each observed system - essentially, the probability, across 
a uniform sample of initial conditions, that the observed pro- 
file was drawn from a simulation with the given starburst mass 
fraction. These are shown, and compared to the fitted extra 
light fraction for our two-component models. In general, the 
fitted extra light fraction corresponds well to the characteristic 

' ' The values shown in Fig ures[lQll4l are based on comparison only to the 
profiles shown, from Kor mendv et aU (2008). In Table[T] the values represent 
the results from all available data sets, including multiple different observa- 
tions of the systems shown here, and so can be slightly different (however the 
differences are generally small). 

The dynamic range of the fits is somewhat difficult to discern in Fig- 
ures [ToIT4] owing to the plotting versus r'/**; we therefore reproduce these 
figures plotting fi versus r in Appendix IaI 
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Fig. 10. — Su rfac e brightness profiles - decomposed into dissipationless and dissipational/extra light components - are shown for core ellipticals in and around 
the Virgo clusteii^l Open circles show the observations, from Kormendy et al. (2008). These are the highest-mass core ellipticals in Virgo (~ 6-8M,). Upper: 
Observed V-band surface brightness profile with our two component best-fit model (solid, dashed, and dotted lines show the total, inner/extra light component, 
and outer/pre-starburst component). The best-fit outer Sersic index, extra light fraction, and rms residuals about the fit are shown. Lower: Colored lines show 
the corresponding surface brightness profiles from the three simulations in our library which coiTespond most closely to the observed system. Dashed line shows 
the profile of the starburst light in the best-matching simulation. The range of outer Sersic indices in the simulations (i.e. across sightlines for these objects) and 
range of starburst mass fractions which match the observed profile are shown, with the rms residuals of the observations about the best-fit simulatiorU-i! Bottom: 
Observed disky/boxy-ness (an) and ellipticity profiles, with the median (solid) and 25—75% range (shaded) corresponding profile from the best-fitti ng si mulations 
above. Note that these are not fitted for in any sense. Figures [T 1414l show the other core ellipticals in the sample, ranked from most to least massivai^. 



Starburst mass fractions in simulations which produce similar 
light profiles. 

6.2. Independent Evidence for the Surviving Dissipational 
Component 

Figure [TSjcompares our estimates of the dissipational mass 
fraction in the observed ellipticals: the directly fitted extra 
light fraction /extra and inferred starburst mass fraction /b 
from the best-fitting simulations. For clarity we restrict to a 
sample where the range in /extra from different sources of pho- 
tometry is less than ^ 20% (we find the same results for the 
entire sample, but with correspondingly larger scatter). We 
compare with the same quantities fit to cusp ellipticals and 
observed gas-rich merger remnants in Paper II and Paper I, 
respectively. Our fitted decomposed extra light fraction re- 
liably traces the inferred starburst mass fraction, with a fac- 
tor ~ 2 scatter similar to that predicted from our simulations 
(Figure |4]i. This is true even for profiles without an obvious 
"break" or feature indicative of a transition to extra light. The 
starburst fraction /sb itself must, in some sense, reflect the cold 
gas mass available in the disks just before the original gas- 
rich merger (and we show in Paper I that this is the case) - in 
this physical sense (insofar as our simulations are reasonable 
qualitative proxies for observed profiles), our fitted /extra and 
inferred /b are a robust reflection of the dissipational content 
of the progenitors. 



We can independently test these decompositions with other 
observations of stellar populations, shown in Figure[T8] Given 
sufficiently detailed observations and stellar population mod- 
els, observers have estimated the mass fraction which formed 
in a more recent, central starburst (as opposed to the more 
extended quiescent star formation history) - completely inde- 
pendent of our analysis here - for s everal of the ob served sys- 
tems j Schweizer 199 6tlTitus et al.|[l9 97; Schweize r & Seitzed 
[1998112007; .Reichardteral.ll200lflMichard .20061) . Unfor- 
tunately, the observational expense required leaves us with 
samples of relatively limited size (there are only a few ob- 
jects with sufficiently deep stellar population data to allow this 
comparison). Nevertheless, comparing our estimated /extra or 
/sb with these estimates for the mass fraction in the secondary 
(newly formed/starburst) stellar populations, we find reason- 
able agreement. They all suggest that our fitted extra light 
component is indeed a good proxy for the mass fraction which 
was involved in the central, merger-driven starburst. We find 
similar good agreement in these cases for both our core ellip- 
ticals here and the cusp ellipticals considered in Paper II, and 
the merger remnants considered in Paper I, together constitut- 
ing a sample of ^ 25 objects. 

Other independent lines of support further suggest that the 
decomposition is recovering a real physical difference in the 
populations. In Figures [T0lfT4l there is often a change in 
isophotal shape and kinematics at the radii we associate with 
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Fig. 15. — Observed surface brightness profiles of a subset of the confirmed core ellipticals from the sample of lLauer e t al.' l'2Q 07a) , with the best-matching 
two component parameterized fit (dashed and dotted lines) and best-fitting simulations (red, orange, and green lines), as in Figure [To] Where multiple sources 
of photometry are available, independen t fits to each are shown. The objects are ranked from brightest to faintest in V-band (as shown): these are among the 
brightest galaxies in the sample; Fi gures fl 6| 1 7 1 continue to fainter luminosities. Profiles are shown over a constant angular scale (top axis; bottom axis shows 
physical radius in kpc). The corresponding (right) panel for each shows the distribution of physical starburst fractions for the simulations which provide a good 
fit to the observed profile (as described in the text), with the fitted (parameterized) extra fight fraction (blue dotted line; one for each source of photometry) 
and observed secondary (recent) starburst components (red dashed, where available). Note that our simulation resolution Hmits do not extend within the central 
~ 30-50pc, and our fits are not intended to describe these radii. 



the dissipational remnant - systems become slightly diskier, 
despite the fact that re-mergers will smooth the isophotal 
shapes and destroy strong transitions in features such as the 
kinematics that might be present in cusp ellipticals at the extra 
light-outer light transition. We show in §[8] that the distribu- 
tion of global kinematic properties and isophotal shapes inde- 
pendently requires re-mergers with a relatively narrow range 
in starburst mass fractions, in agreement with our estimates 
from fitting the surface brightness profiles. 



Moreover, we show in Paper 1 that the extra light is asso- 
ciated with, and indeed drives, stellar population gradients in 
the remnant elliptical. In particular, we argue that metallicity 
gradients are most robust (age gradients are very weak once 
the system is old, and gradients in a-enhancement are much 
more sensitive to initial conditions and the exact details of e.g. 
the starburst timing relative to the timescale for star formation 
in progenitor disks). We show in §|5]that these are expected to 
survive re-mergers relatively intact, and therefore should pro- 
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vide independent indicators of the dissipational versus dissi- 
pationless decomposition. 

There are a small number of systems in our observed core 
sample for which detailed stellar population gradients have 
been measured in ^anchez-Blazquez et al. (2007) or two- 
dimensional stellar pop ulation maps have been constructed in 
iMcDermid et alj (120061) . In almost every case, we see what is 
predicted: a strong rise in metallicity around the radii where 
our fits infer the transition between the dissipational and dis- 
sipationless components, with the metallicity gradient flatten- 
ing within the radius where the extra light dominates the pro- 
file. Figure [19] shows a few illustrative examples, compar- 
ing the observed metallicity gradients to our two-component 
surface brightness decompositions. If we attempt to quantify 



where the metallicity profile flattens (expected inside the ra- 
dius where the extra light dominates, since the predicted gra- 
dient is largely driven by the transition from an outer, less 
metal-rich population into the more metal-rich dissipational 
population) by fitting it to an external power law with a break 
at some radius Rt, and compare this radius Rb to the radius 
at which the fitted inner/starburst component dominates the 
profile (the outer radius where the fitted inner component sur- 
face brightness becomes larger than that of the outer compo- 
nent, or similarly the effective radius of the extra light com- 
ponent), we find a remarkably good correlation. In every 
case, the age gradients (while weak and less uniform, as pre- 
dicted) also confirm that the central component we identify 
is younger than the outer profile, as predicted (the difference 
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is usually small, 1 Gyr, as expected in merger models with 
relatively early formation - if the central components were in- 
stead formed by e.g. dissipation of gas from stellar mass loss 
(ICiotti & Ostriker 2007), we might expect the central compo- 
nents to be much younger). More data are needed to test these 
decompositions and predictions over the entire dynamic range 
of interest, but our preliminary comparisons all lend some 
confidence that we are, in an average sense at least, recovering 
the appropriate physical separation between inner radii domi- 
nated by the remnants of extra light systems originally formed 
in dissipational starbursts and outer light profiles formed by 
violent relaxation acting on older, initially low-phase space 
density (disk) stars. 



7. PROPERTIES AND SCALING LAWS OF DISSIPATIONAL 
COMPONENTS 

7.1. Amount of Dissipation Versus Mass 

The preceding analysis indicates that re-mergers do indeed 
have a "dissipational" component (the surviving remnant of 
the original starburst that formed the progenitors' extra light), 
and that its mass can be estimated from their surface bright- 
ness profiles. Now, we examine how this scales with other 
properties, specifically in reference to how the extra light'^ 

" In what follows, we use the phrase "extra light" for the core ellipticals in 
the physically pertinent sense, to refer to the component that traces or reflects 
the oiiginal starburst mass, i.e. the centrally concentrated contribution to the 
light profile that is the remnant of the original gas-rich merger remnant extra 



Dissipational Signatures in Core Ellipticals 



21 



1.0 



0.1 




1.0 



0.1 




1.0 



js 0.1 




1.0 



J 0.1 



Independent SED-Based Estimates 
□ 



0.1 1.0 
starbursi (best-fit simulatioRs) 




(stellar populations) 



«- 1.0 



tZ 0.1 




«- 1.0 



tn 0.1 



Independent SED-Based Estimates 



•tfr Kormendy et al. 
O Lauer et al. 
IT Rothberg & Joseph 




0.1 1.0 

starbu,si(stellar populations) 



^^^^ Tactio n in the best-fitting simu lations 

(/st; as in Figures [lOllSK We show the results from the core elliptical sample studied here (bottom; red stars are systems from lKormendv et alj (2008"), violet 
circles from Lauer et al. (2()07a)) th e obse rved gas-rich merger remnant and cusp elliptical samples in Paper I and Paper II (top; cyan squares are the merger 
remnants from Rothberg & JosepH 120041) ). The two estimates agree well, with a factor ~ 2 — 3 scatter in /cxtra(/sb) (similar to what we expect from our 
simulations; see Figure|4). Open point with error bars shows the typical range in /extra and /sb across different sources of photometry and in different wave bands. 
Center: Fitted /extra versus independent observational estimates of the mass fraction formed in a more recent starburst/star formation event, from two-component 
stellar population model fits to the observed SEDs iSchweizei 1996; Titus et al. 1997; Schweizer & Seitzer 1998, 2007; Reichardt et al. 2001; Michard 20061). 
Right: Same, but comparing the stellar population estimates to /sb from the best-fitting simulations. More observations are needed to test our estimates, but the 
stellar population data independently suggest that our decompositions are reasonable. 



V 



NGC 1700 




NGC 3379 




NGC 4472 



1.0 - 



NGC 1600 



0.01 0.1 



0.01 0.10 1.00 
r [kpc] 



0.1 1.0 
r [kpc] 



0.1 1.0 10.0 

r [kpc] 



0.1 
0.1 





/ 

y 

/ 




|- / 


1 1 


✓ 




✓ 

/ 


1 jsf 


1 1 
✓ 









1.0 

R (metaiiicity profile break) [l<pc] 



Fig. 19. — Independent tests of our two-component decompositions with observed stellar population gradients in core ellipticals. Left: Comparison of the 
surface brightness profile (ton; points) and two-component decomposition (dashed and dotted lines show the fitted inner & outer components for each source 
of photometry, as Figure [Tsl with the observed metalHcity profile (bottom; as discussed in text this is the most robust indicator of extra light) for the cases 
observed in Sanchez-Blazguez et alj 120071) . The comparison agrees with our predictions in Paper II and §[5] the metaiiicity gradients are strong at the radii of 
the transition from the outer (dissipationless) stellar populations to the inner (dissipational, metal-enriched) population, and begin to flatten within the radius of 
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content was seen to scale with galaxy properties in cusp ellip- 
ticals in Paper II. 

Figure |20] plots the inferred starburst mass fraction f^\, for 
the observed systems as a function of stellar mass. In the same 
manner that we have defined a best-fit /sb from the best-fit 
simulations, we can also define a best-fit "initial" gas fraction 
(roughly the gas fraction ~ I Gyr before the final merger), and 
show this as well. We emphasize though that this is a much 
less robust quantity (for example, if we place the systems on 
a longer period orbit with the same initial gas fraction, then 
by the time they merge, more gas will have been consumed in 
the disks and the results will look like a shorter period orbit, 
lower "initial" gas fraction simulation). In either case there 
is a clear trend of increasing dissipation (increasing fractional 
mass required in a dissipational starburst component) at lower 
masses. We compare this with the behavior found for ob- 
served cusp ellipticals in Paper II. The significance of the cor- 
relation is unambiguous (> Sfi, in each of the cusp and core 
samples). The two trends are statistically indistinguishable: 
fitting them as a function of stellar mass with the form 

(/sta.Wt)~ J , (3) 

we obtain (Mo, a) = (IQ'^-^^^-^Mq, 0.43 ±0.04) for the cusp 
sample and (10** '*='=°-^Mq, 0.35 ±0.04) for the core sample, 
with roughly a constant factor 2 intrinsic scatter at each 
mass in both cases. It is worth noting that if identical systems 
begin on the cusp correlation and re-merge, they will conserve 
/starburst but doublc M, . To the extent that there is not a large 
offset between the two correlations, then, this suggests that 
at least the typical core ellipticals have not undergone a large 
number of re-mergers. A single dry re-merger, especially if 
it is of a more typical ^1:3 mass ratio, is however consistent 
with what we see, within the scatter 

In Paper II we noted that the inferred gas fractions corre- 
spond well to observed disk galaxy gas fractions as a func- 
tion of mass, from z = 0-2. We show these observations as 
well, and note that they are sufficient to explain the trend and 
scatter in the extra light fractions of both cusp and core el- 
lipticals. In other words, the distribution in progenitor gas 
fractions implied by the elliptical surface brightness profiles 
is, as a function of mass, exactly what would be predicted if 
one assumes that the original progenitors were spirals, and 
that most of the systems were first formed by a major merger 
somewhere between a redshift of z ^ 0-3. In fact, this is 
exactly what is inferred for the formation times of the ellip- 
ticals in our sample and those of similar mass (< a few L*) 
from both observations of the early-typ e or red galaxy mass 
funct i ons dBundv et a l. 2006; Borch et al .|2006t[Fontana et alj 
llOOe"; 'Hopkins et al.' 2006b, 2007a, 2008c), from direct stellar 
population synthesis studies (Trager et al. 2000; Thoma s^t alJ 
I2OO5I: iGallazzi et al.l 120051) . and by association of elliptical 
galaxy formation with the triggerin g of quasar activity (e.g. 
iHopkins et al.ll2006ai 12007112^83) . 

There may be an interesting suggestion that the extremely 
massive, BCG tail of the core population (for which the vari- 
ous caveats described in §|4]apply) might asymptote to ~ con- 
stant values of /sb ^ 0.05-0.1 atM, > lO'^M*. This would 
be expected if these systems all "began" at lower masses 
(where this gas fraction is typical of disks) and moved up by a 
large number of dry re-mergers. Unfortunately, the data here 
are limited, but we stress that these systems do not signifi- 
cantly alter our comparison, and are not the systems we in- 
tend our modeling (applicable to those core ellipticals with a 



moderate re-merger history) to describe. 

The inference of /sb requires some comparison with our 
simulations, but /extra is a purely empirical quantity. We there- 
fore repeat this exercise with the empirically fitted extra light 
component /extra, and show the results in Figure|2T| The trend 
seen in /extra(A^*) is completely consistent with that in /b, but 
with a scatter larger by a factor ^ 2 (exactly what we expect, 
based on the predicted and observed scatter in /extra(/sb))- 
Moreover the trend we find in the core ellipticals is consis- 
tent with that in the cusp elliptical sample from Paper II. Con- 
sidering just the data in Figure [21] even given its increased 
scatter, the trend of decreasing extra light fraction with mass 
is significant at > 6 ct (> 8 cr, if we include both the cusp and 
core ellipticals). We have experimented with alternative, non- 
parametric (albeit less accurate) estimators based on e.g. the 
concentration indices or stellar populations of our simulations 
and observed systems, and obtain a similar answer In short, 
even without reference to our simulations, however we derive 
an estimate of the dissipational component, it is difficult to 
escape the conclusion that it is more prominent in lower-mass 
ellipticals, even in core ellipticals that may have experienced 
re-mergers. The mass range covered by the two samples is 
somewhat different, but where there is overlap they are simi- 
lar, and at the extremes the two populations appear to be con- 
tinuous extensions of the same trends. 

7.2. Sizes of the Extra Light/Dissipational Region 

The inferred dissipational components in simulations, cusp 
ellipticals, and core ellipticals also appear to follow a similar 
size-mass relation (especially if we restrict to systems with 
similar /sb), shown in Figure |22] In Paper I we show that the 
size-mass relation is driven by the condition that the gas col- 
lapsing into the central regions in the final starburst becomes 
self-gravitating: (GMextia/^extra) (GM*//?,). In Figure |22] 
we show that the observed core ellipticals also follow this con- 
straint, with small scatter and dynamic range similar to that in 
our simulated mergers. This is expected even if these sys- 
tems have experienced substantial re-merging, provided that 
re-mergers puff up the system uniformly. 

To test, this. Figure |23] plots the distribution in the ratio 
of the effective velocity dispersion/self-gravity of the inner 
component to that of the outer component, for the cusp and 
core elliptical samples. Interestingly, there is small offset, in 
the sense that the effective dispersions of the core ellipticals 
are slightly higher than those of cusp ellipticals. This is true 
even if we restrict the analysis to a narrow mass bin - i.e. 
the extra light in core ellipticals is (on average) slightly more 
concentrated relative to the outer light component, compared 
to cusp ellipticals. 

It is unlikely that this is caused by scouring or whatever 
process might directly give rise to the core on < 50pc scales, 
since the effect is seen even in the high-mass objects where the 
effective radius of the extra light component is > 1 kpc. How- 
ever, it is well established in numerical simulations that, when 
merging coUisionless systems, angular momentum and orbital 
energy are preferentially transferred to less bound material, 
expanding their orbits while allowing the central reg ions to re- 
main tightly bound (see, e.g. i Hernquistll 1 992l 1 1 993i\ In prac- 
tice, most of this less bound material is part of the outer halo 
- so the gradient in terms of energy input and expansion of 
material across the baryonic galaxy is weak (i.e. the profiles 
are stretched by a roughly constant factor in Figure[Tll. Nev- 
ertheless, a small effect of this nature would be expected in 
coUisionless re-mergers. The observed offset is in fact weak. 
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Fig. 20. — Inferred gas content (dissipational/starburst fraction) of elliptical-producing mergers as a function of stellar mass. Initial gas fraction (top) and 
physical final starburst mass fraction (bottom) corresponding to the best-fit simulations to each observed system in the samples of ILauer et a l.' K007aO (circles) 
and Kormendy et al. (2008) (stars) are shown, with the typical 25 — 75% allowed range (eiTor bar). We show results separately for cusp ellipticals (left), core 
ellipticals (center), and both together (right). Dashed (solid) line shows the fit to the data (Equation [3) in cusp (core) ellipticals. Colored points with error 
bars indicate the mean (and ± 1 cr range in) disk gas fractions at the same stellar mass, at = I Bell & de Jong 2001 ; Kannappan 2004; McGaugh 2005, blue 
diamonds, squares, and circles, respectively), z = 1 I Sh apley et al. 2005. green squares), and z = 2 lErb et al. 2006, orange circles). There is a clear trend of 
increasing dissipation required to explain elliptical profiles at lower masses (significant at > Scr in either core or cusp subsamples separately), in good agreement 
with the observed trend in progenitor disk gas fractions over the redshift range where cusp elli pticals are f ormed, and with what is invoked to explain the observed 
densities and fundamental plane correlations of ellipticals (e.g. Kormendy 1989; Hernguis t et alJI 19931) . The best-fit trends in cusp and core populations are 
statistically identical; i.e. the dissipational/extra light component is preserved in re-mergers. 



a factor ~ 1.5 in /^extra/pouter (in the most extreme case of 
a single major re-merger where the inner profile was not af- 
fected by the merger and the outer profile doubled in radius, 
we would expect a factor of two offset), and more data would 
be needed to demonstrate that it is robust. However, if real, it 
provides further support to the idea that these have re-merged, 
but not destroyed their extra light content. 

7.3. Dissipation and the Sizes of Ellipticals 

We demonstrated in Paper II that the effective radii of cusp 
ellipticals scales with their extra light content - systems with 
more dissipation are smaller, since more of the light comes 
from the compact starburst remnant. Figure |24] extends this 
to the core elliptic als. Specifically, we determine {Rg{M^,)) 
from the sample of lShen et al.l (l2003l) . and take the ratio of the 
half mass radius of each system (determined directly from the 
light profile, or from the fits, it does not change the compari- 
son) to that value. Our mass bins are small enough, however, 
that this makes little difference compared to just e.g. consid- 
ering Re in a given bin. Again, we see a similar behavior in 
the cusp and core ellipticals: at a given stellar mass, systems 
with larger extra light have systematically smaller (they 
also have slightly larger velocity dispersion cr, although the 
scatter is larger there in both simulations and observations). 
We showed that this could be predicted analytically, follow- 
ing |Covingtonetal] ( l2008h given the impulse approximation 



to estimate the energy loss in the gaseous component, fol- 
lowed by collapse in a self-gravitating starburst. For typical 
conditions, this reduces to the simple approximation 

Reifsb = 0) 



R. 



(4) 



l + (/sb//o)' 

where /o « 0.25-0.30 and Reifsh = 0) is the radius expected 
for a gas-free remnant. In each case, the simulations and ob- 
served systems occupy a similar locus to that predicted in this 
simple dissipational model. We can also construct this plot 
with the starburst mass fraction /sb of the best-fitting simu- 
lation as the independent variable, and find an even tighter 
correlation of the same nature. Core elliptical sizes depend 
on the relative strength of dissipation (evident in their profile 
shapes) in the same manner as in cusp ellipticals. 

This is expected: even if the systems with cores have ex- 
panded via re-mergers, they should (so long as there is not 
a wide range in number of re-mergers or systematic depen- 
dence of the number of re-mergers on starburst fraction at 
fixed mass) grow uniformly, and preserve these trends. We 
might expect some normalization offset: if the mean size- 
mass relation after a gas-rich merger is a power law R,, oc M", 
and two such systems with mass ratio / (where / is the mass 
ratio of the secondary to the primary) have a parabolic dry 
merger (and preserve profile shape), then the remnant will in- 
crease in radius by a factor (1 -l-/)^/(l +/^~") relative to the 
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Fig. 21. — As Figure l20l but showing our empirically fitted /extra as a 
function of stellar mass. The trend of increasi ng d issipation at lower masses 
is still clear and is consistent with that in Figure l20l but with an expected extra 
factor ~ 2 — 3 scatter from the scatte r in our purely empirical estimator. Solid 
line shows the best-fit from Figure |20| dotted lines sh ow the ± 1 cr scatter 
expected based on the scatter in /c.xtra(/sb) (see Figure [Ts). Fitting to just 
these data yield trends statistically consistent with the solid line, ruling out 
no dependence of dissipation on mass at > 8 cr confidence. 

primary. However, it has also increased in mass, so compared 
to ellipticals of the same (final) mass, its relative increase in 
size is only (1 + /)^~°/( l +/^~"). For observationally sug- 
gested values a = 0.56 (Shen et al.ll2003 ). this yields only a 
30% (~ 0.1 dex) relative size increase for mass ratios of 1:3 
through 1:1. This is easily dwarfed by the effects of dissipa- 
tion seen in Figure|24l which can change the sizes of systems 
by nearly an order of magnitude. That there is not a sub- 
stantial normalization offset between the trends for cusp and 
core ellipticals in Figure |24] therefore suggests that the num- 
ber of major re-mergers for typical core ellipticals has been 
relatively moderate, but it is not a strong constraint. Our simu- 
lations and the observations plotted suggest that the dominant 
factor regulating the size of an elliptical is the amount of dis- 
sipation, even for systems which have undergone re-mergers. 

8. DISSIPATION AND RE-MERGER KINEMATICS/SHAPES 

In Paper II we showed that the existence and mass fraction 
in a dissipational component is related to the kinematics and 
shapes of cusp ellipticals. In particular, we confirmed what 
has been recognized in a growi ng number of simulations (see 
ej. Naabetal. 2006at ICoxet al. 2006b; Jesseit et al. 2007; 
iBurkert et al.l l2008h . that disk-disk mergers without signif- 
icant dissipational components (whether this component is 
created in the merger, as a gas-rich merger, or "inherited" 



from the progenitor disk bulges - which formed it dissipa- 
tionally in some earlier encounter/dissipational event) do not 
resemble any observed spheroids: they are slowly rotating and 
boxy, but highly elliptical. However, systems with the appro- 
priate dissipational content for their stellar mass (as in e.g. 
Figure |20] l agree well with observed distributions in cusp el- 
lipticals in all of these properties (and the same is true for 
more sophisticated kinematic and shape measures, includ- 
ing their orbital isotropy, rotation, triaxiality, kinematic mis- 
alignments, and hig her-order velocity moments /13 and I14; 
ICox et al. 2006b; Jesseit et al .| [200 7). This is related to two 
aspects of the gas content in the mergers. First, the dissi- 
pational stellar component itself, by creating a dense central 
stellar mass concentration, acts for stars at almost as a cen- 
tral point mass with large stellar mass fraction ^ 10-20% - 
directly making the galaxy potential more spherical (less el- 
liptical) an d reducing the portion of phase space available to 
box orbits CBarnes & HernquistI 199 111 19961) . Second, if more 
gas is available to form the starburst, more is likely to survive 
the merger This gas will quickly re-form a disk (albeit often 
a very small, nuclear disk with relatively small mass fractions 
of order a couple percent, more often designated a nuclear 
subcomponent, if it can be recognized in observations), giv- 
ing rise to more rapid rotation and diskier isophotal shapes. 

Given the connection between the dissipational component 
and the shapes/kinematics of the disk-disk merger remnant, 
it is worth examining how this compares after dissipation- 
less spheroid-spheroid re-mergers (and how this scales with 
the original starburst/dissipational mass fraction). Figure |25] 
shows this comparison for representative simulations with dif- 
ferent gas fractions. We consider the rotation (V /a, where V 
is the maximum rotational velocity from the projected stel- 
lar rotation curve and a the central dispersion), ellipticity e, 
and isophotal shape (boxy/diskiness, quantified as 04/0 in 
standard fashion) for a 1:1 identical disk-disk merger and 
spheroid-spheroid re-merger of that remnant 

In the case where the initial disk-disk merger is gas-poor 
(/starburst ~ 0.01), wc rccovcr the results in the papers de- 
scribed above: gas-poor disk-disk mergers yield remnants 
that are slow rotators with large ellipticity, without a corre- 
lation between e and 04, and extending to very boxy isopho- 
tal shapes, in stark disagreement with the properties of either 
cusp or core ellipticals (given similar resolution and seeing). 
This is relatively unchanged in re-mergers, and remains so re- 
gardless of the number of "generations" of re-mergers (we 
have experimented with series of several such re-mergers): 
without dissipation to build a central mass concentration (in- 
crease the central phase-space density), dissipationless pro- 
cesses cannot produce such a central density as required to 

Specifically we show the rotation (V/c) or (V /a)* measured within 
Re, and boxy/disky-ness lOGa^/a. and ellipticity e measured for the half- 
mass projected (R^) isophotal contour, for each of 100 lines-of-sight to 
the remnant uniformly sampling the unit sphere (i.e. representing the dis- 
tribution across random viewin g angle s ). T he details of the fitting for our 
simulations are described in Cox et al.' f2006b), and for the observations in 
[Bender et al. 1 1987); Bender 1 1988); Bender et al. ( 1992); Faber et al. 1 1989); 
[Simien & Prugniel ( 1997a, 2002). We define rotation in the standard manner, 
relative to that of an oblate isotropic rotator, with the parameter (V/cr)*, de- 
fined using the m aximum circular velocity Vc, velocity dispersion within i?e, 
and ellipticity as jKormendvl 1 982ft 

(V/aT =(V/a)l^el(\-e). (5) 

We exclude the coplanar merger simulations from our comparisons here; 
those simulations are idealized perfectly coplanar prograde orbits, and as such 
produce pathological orbit and phase structure (we do, however, include some 
more representative orbits below that are not far from coplanar). 
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Fig. 22. — Top: Effecti ve r adius of the extra light component (not equivalent to the radius where it breaks from the outer Sersic fit) as a function of extra Ught 
mass (points as in Figure [Tsl . Simulations (left) and observed cusp (middle) and core (right) elhpticals are similar, especially if we restrict to simulations with 
initial gas fractions ~ 0.2 — 0.4. Resolution limits (see Paper II) prevent us from simulating systems with Scxtra ^ lOOpc, but this is only important for the few 
very lowest-mass ellipticals (L <0.01L,). Filled diamond is typical sightline-to-sightline variance in the simulations, open diamond the source-to-source (or 
band-to-band) scatter in observed profile fits. Bottom: Effective velocity dispersion of the extra light component vs. that for the whole galaxy. Solid line shows 
(G Msxim / Rcxtm) = (GMt/Re) - the extra light collapses to the point where it is self-gravitating. The observed systems follow a trend which agrees well with the 
simulations. 
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Fig. 23. — Cumulative (top) and differential (bottom) distribution of the 
ratio of effective velocity dispersion of the inner light component to that of 
the outer light component ((GMcxtia/Scxtra)'^"/(GM,/Rf )'/2), for cusp el- 
lipticals (blue) and core ellipticals (red). Dotted line shows the median for 
each. Typical values are similar, but slightly higher for core ellipticals (the 
difference in mean/median is significant at > 3 cr), suggesting that the inner 
component is relatively less puffed up by re-mergers than the outer compo- 
nent. 



make the potential sufficiently round and eliminate most of 
the phase space available to the most extreme boxy/triaxial or- 
bital configurations (which are not observed in real systems). 
Regardless of re-merging, initially gas-free disk-disk mergers 
cannot ultimately reproduce the kinematics and shapes of core 
ellipticals. 

If instead the initial disk-disk merger is reasonably gas-rich 
(/starburst ~ 0.1), in the range that we have inferred from fit- 
ting observed surface brightness profiles of ^ 1 - 5 ellip- 
ticals (and in agreement with the observed gas fractions of 
disks in the relevant mass and redshift range), then the rem- 
nants are more round (owing to dissipation building a central 
mass concentration), with kinematic subcomponents that give 
rise to the rotation and disky isophotes in good agreement 
with the observed cusp population (see e.g. iKrajnovic et all 
I2OO8I who show such subcomponents are ubiquitous in cusp 
ellipticals and drive much of this rotation). Re-mergers of 
these remnants destroy these kinematic subsystems (they are 
small disks/subsystems, and as such are violently relaxed in 
re-mergers), but retain the central mass concentrations (as we 
have argued they must), and this leads to the formation of 
sufficiently round, but boxy and slowly-rotating systems with 
the e-fl4 correlations appropriate for core ellipticals (com- 
pare the observa tions in e.g. Bender 1988; Bender et al. 19921; 
lEmsellem et alJ i2004. 2007). The re-merger remnants of cusp 
ellipticals with the starburst mass fractions we have indepen- 
dently inferred from profile fitting appear to agree well with 
observed core elliptical kinematics and isophotal shapes. 

If the initial disk-disk merger is too gas-rich (/starbunst ~ 
0.4-0.6), larger than that inferred for any of the observed core 
ellipticals in our sample, this presents an interesting case as 
well. These extremely gas-rich situations may apply to some 
of the most low-mass cusp ellipticals, but are clearly more 
gas-rich than expected at the masses of typical co re ellipti- 
cals (and disk progenitors), even at high redshifts jErb et alj 
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Fig. 24. — Effective radius Re relative to the median value for all ellipticals of the same stellar mass, as a function of our fitted extra light fractions (the 
empirical tracer of the dissipational/starburst mass fraction). We compare simulated merger remnants (top) with observed cusp ellipticals and gas-rich merger 
remnants (middle) and observed core ellipticals (bottom), with points as in Figuresl4l& ll8l We show this in three bins of stellar mass (relative to M, lO" Mq , 
or My = —21). Solid (dotted) lines show the mean (ilcr) correlation, following the analytic solution for dissipational mergers and fits to our simulation in 
ICovington et alj 4200 8). Simulations and observations exhibit the same behavior: systems with smaller Re at fixed mass have systematically higher extra hght 
fractions (> 6a significance in the observations). This implies that, at fixed mass, systems are driven along the fundamental plane by the relative amount of 
dissipation involved in their formation. This behavior is true regardless of cusp/core status, suggesting that (most) core ellipticals have not had a large number of 
re-mergers. 



|2006^. In these cases, we find that the initial disk-disk merger 
remnant is even more round, rapidly rotating, and disky. How- 
ever, with such a large dissipational component, the kinematic 
subsystems are often quite massive - in many cases they are 
more appropriately d escribed as large surviving disks (see e.g. 
iHopkins et al.' 2009). As a consequence, it becomes difficult 
to simply disrupt them in re-mergers. Moreover, conserva- 
tion of the specific angular momentum in such large compo- 
nents makes it very difficult to achieve a slow rotator even 
if the subsystem is disrupted (as opposed to the case above, 
where the relevant subcomponents represent a small fraction 
of the galaxy mass); this has been seen when e.g. mock el- 
lipticals with uniform rotation (as opposed to more realistic 
rotation associated with kinematic subcomponents) are re- 
merged (iNaab et al.ll2006bt iBurkertet aLll2()08t iJesseit et all 
l20(MlCox etal.l2009 ()~ As a consequence, the re-merger rem- 



nant looks like an oblate isotropic rapid rotator, not a slow 
rotator, and exhibits correlations between e.g. diskiness and 
ellipticity (and higher-order kinematic moments) that do not 
resemble either the cusp or core population. 

Some of the observed differences between cusp and core 
distributions owe to the fact that core ellipticals (because 
they predominate at high masses) tend towards lower /sb 
and cusp ellipticals (dominant at low masses) towards higher 
/sb. However, at the same mass, the two populations evi- 
dence similar /b (Figure |20| |. but show differences in their 
kinematics. Figure |25] demonstrates that at fixed /sb the 
kinematics also depend significantly - and in a different 
manner - on merger history. As a consequence, the well- 
established observationa l diffe r ences between the se popula- 
tions (e.g. iBender et all 119891: iFaber et al.1 Il997h are non- 
trivial - there is no population with a large number of dis- 



Dissipational Signatures in Core Ellipticals 



27 



> 



1.2 
1.0 

0.8 
0.6 
0.4 

0.2 

0.0 

0. 



Gas-Poor (f^i, < 0.05) 

o Dissipational 
(Disk-Disk) Merger^ 

□ Dissipationless 
(Spheroid-Spheroid) Re-Merger 




ni 

o 
o 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Ellipticity e 




o 




eo 03— 

' a 

« Observe 


BB o g 

d 



.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Ellipticity e 



-1 1 
1 00 aja 





> 




03 

O 
O 



.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Ellipticity e 


















fgt, ~ Observed 





.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Ellipticity e 



-1 1 
1 00 aVa 



> 




to 

o 
o 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Ellipticity e 








» Observe 


d 



.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Ellipticity e 



-1 1 
1 00 aVa 



Fig. 25. — Rotation (V/cr), ellipticity e, and isophotal boxy/diskiness (a4/a) for a 1:1 identical disk-disk merger and spheroid-spheroid re-merger of that 
disk-disk remnant. Top: Case where the initial disk-disk merger is gas-poor (/staibmst ~ 0.01). Middle: Case where the initial disk-disk merger is gas-rich 
(/stiuburst ~ 0. 1, typical of ~ 1 — 5L, ellipticals). Bottom: Case where the initial disk-disk merger is extremely gas-rich (/staiburst ~ 0.4 — 0.6, larger than we infer 
for any of the observed core ellipticals in our sample). For each, we show a single merger remnant, as viewed from ^ 100 independent sightUnes. Gas-poor 
disk-disk mergers yield remnants that ai'e slow rotators with large ellipticity, without a correlation between e and a^, and extending to veiy boxy isophotal shapes. 
This is relatively unchanged in their re-mergers (simple energetic concerns mean they cannot get significantly "more round," and the extension to "too boxy" 
isophotes and lack oit—a^ correlation is enhanced). Gas-rich disk-disk mergers ai'e more round (owing to dissipation building a central mass concentration), with 
kinematic subcomponents that give rise to the rotation and disky isophotes observed in cusp elHpticals (and their correlations). Their re-mergers destroy these 
subsystems, but retain the central mass concentrations - this enables the formation of veiy round, but boxy and slowly-rotating systems with the e-a^ coiTelations 
appropriate for core ellipticals. But systems which are significantly more gas-rich than we infer remain disky and rapidly rotating even after re-mergers, in conflict 
with the observed core population. 



sipationless re-mergers (at any /sb) that matches the observed 
kinematic/shape distributions of cusp ellipticals, and no pop- 
ulation of disk-disk mergers (at any /sb) that matches the ob- 
served distributions of core ellipticals'^. 

Figures|26]&|27]show how the rotation and isophotal shapes 
of our simulated disk-disk and re-merger remnants compare 
with observed cusp and core elliptical populations. We specif- 
ically consider simulations with similar dissipational fractions 
to the observations (we only compare simulations that have 
values of /sb in broad agreement with the locus we observa- 
tionally infer from fitting the observed systems in Figure l20l 

" As discussed below, this statement is not true for any single object, but 
applies to the populations as a whole. 



equivalently but more strictly, we can consider only the set of 
simulations each of which is a best-fit match to one of the ob- 
served surface brightness profiles, as in Figure [TOb. Despite 
these simulations being constrained to match the observed 
systems only in their averaged surface brightness profile, they 
occupy a similar locus to the observations in each sense in 
Figures |26]&|27] 

Figure |27] illustrates that the distribution in e, 04 /a and 
(y/fj)* for core ellipticals constrains not only their merger 
history, but also the dissipational content of the original 
gas-rich mergers, (even if these systems are the product of 
spheroid-spheroid re-mergers). This constraint agrees well 
with our estimates from their surface brightness profiles: sys- 
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Fig. 26. — As Figure l25l but comparing our library of simulations to observed systems. Top: Rotation {V/a), ellipticity e. and isophotal boxy/diskiness (04/0) 
for disk-disk merger remnants, with dissipational/starburst fractions appropriate for their mass (matched to that infen'ed from the observed systems; roughly 
/starbuist ~ 0. 1 for most of the near ~ L, sample). Each point is the median value for a given simulation, across ~ 100 sighthnes. Below, we compare observed 
cusp ellipticals and gas-rich merger remnants. Bottom: Same, for spheroid-spheroid re-mergers of the original disk-disk merger re mnan ts, and observed core 
ellipticals. The simulated and observed locii agree for both populations (despite being distinct from one another). Comparing to Figure l25] it is clear that mergers 
and re-mergers of systems without sufficient dissipation (less than that inferred from our surface brightness profile fitting) do not agree with the observed cusp 
and core populations - significantly different dissipational fractions from our estimates are ruled out by the kinematic data. 
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Fig. 27. — Histograms show the distribution in each property from Figure l26l (ellipticitv e. isphotal shape boxy/diskiness 04/0, and rotation (V /cr)*), for 
observed systems and simulations. Top: Observed cusp ellipticals (black line) compared to simulated disk-disk merger remnants which have similar moderate 
extra light fractions (blue shaded), and gas-poor (nearly dissipationless) disk-disk mergers (/si, < 0.1; red dotted). The simulation distributions uniformly 
sample each simulation in solid angle over ~ 100 lines-of-sight - sightline-to-sightline and object-to-object differences contribute comparably to the scatter. 
Bottom: Same, for observed core ellipticals (black line) compared to simulated re-mergers of the disk-disk merger remnants with appropriate extra light fractions 
(dissipational fractions in the original disk-disk merger; blue shaded) and with much lower (red dotted) or higher (green dotted) dissipational content. The 
simulations, when they have the appropriate extra light/dissipational content, agree with the observed distributions. Their is some "memory" of the amount of 
dissipation retained in the kinematics; systems without sulScient extra light are too elliptical and have a much broader 04 distribution extending to extremely 
boxy isophotal shapes; systems with too much extra light remain too disky and rapidly rotating even after re-mergers. 
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tems which originally had too much or too little gas, relative 
to that we infer from the surface brightness profiles, fail to re- 
produce the observed shapes and kinematics (with or without 
re-mergers). 

Unlike the case of surface brightness profiles, however, 
such a constraint can only be derived from a statistical popu- 
lation of ellipticals; as in Figure |25] the sightline-to-sightline 
scatter (even the object-to-object scatter in median properties 
seen in Figure|26]l is sufficiently large that the isophotal shape 
and/or rotation of any individual elliptical, observed from 
only the single available viewing angle, does not present a 
strong constraint on the gas-richness of the original spheroid- 
forming merger (especially after allowing for the dependence 
on merger history in a re-merged population). However, given 
the distribution observed in a statistical population, it is pos- 
sible to make a strong statement of independent support for 
typical dissipational fractions similar to that inferred from the 
individual surface brightness profile fits. 

Finally, we note that one population that our simulations 
may have difficulty reproducing are the most extreme slow- 
rotators, with (V /a)* < 0. 1 . This may be a by-product of ob- 
servational uncertainties or differences between our simula- 
tion versus observational fitting methodologies. At this level, 
both simulations and observations are sensitive to resolution 
and to the outer radii within which velocities are defined (with 
velocities, eve n in slow rotators, risi ng significantly at larger 
radii; see e.g. lEmsellem et aT]l2007l) : many of the observed 
objects in this category have V/a measured within relatively 
small radii < 0.6 Re- We have broadly matched the simulation 
resolution and radii where kinematic quantities are defined 
(= Re) to most of the observations shown, but a more rigor- 
ous comparison should carefully match the sampling radii for 
given sub-classes of objects and consider full kinematic pro- 
files, rather than the integrated quantities here. 

The difference between our simulations and these slowest 
rotators may also reflect a real physical difference: probably 
our simplific ation of realistic merger his t ories. As demon- 
strated in e.g. lBovlan-Kolchinet"a!T ( l2005L[2006h . the rotation 
properties of spheroid- spheroid merger remnants can be sig- 
nificantly dependent on orbital parameters. We do not find as 
strong a dependence as those authors (namely because their 
model ellipticals did not include a central mass concentration 
as do the systems here; the presence of this component leads 
to more regular tube and loop orbits as opposed to the long 
box orbits that support the most extreme ellipticity/rotation 
configurations seen in that paper); however, we still expect 
that a systematic bias towards certain orbital parameters (for 
example, a preference for particularly radial orbits in the for- 
mation of the most massive ellipticals, if they accrete galaxies 
along filaments) can change (V/a)* in a corresponding sys- 
tematic manner. 

Also, the systems with such low rotation are typically 
the most massi ve systems, with early s t ellar formation 
times ( see e.g. iBender etalJ 1199^ 119931: iFerrarese et alJ 



momentum and such littl e rota tion (iGallagher & Ostrikei) 
1972t[O striker & Tremain dll 1 975t iHa usman & OstrikerUaTT 



Weil & Hernauist.l994,1996HNaab et al..,2007,) . The forma 



1994 'Kuntschner eFal l2006t lEmsellem et al J 120071: 
Cappellari et al. 200^ Kormendv et al. I200a) : these are 
most likely to have undergone a mix of major and minor 
mergers, with e.g. five 1:5 mergers becoming a more likely 
channel for mass doubling in dry mergers than a single 1 : 1 
merger for the most massive systems L ^ . Although this 
will not change our analysis of the extra light components in 
the surface brightness profiles (provided these are as robust to 
minor mergers as they are to major mergers), it has long been 
proposed as a means to yield systems with low net angular 



tion of these objects, therefore, requires further study, with 
more cosmologically representative accretion and merger 
histories. 

9. OUTER SERSIC INDICES OF ELLIPTICALS WITH CORES 

Given the issues surrounding profile fitting described in §|4l 
we have to be somewhat careful when we attempt to describe 
the typical outer Sersic indices of core ellipticals. We can 
(and do) obtain different answers depending on how we ap- 
proach our fits, even for data sets of excellent quality and 
large dynamic range. However, we still find that the distri- 
bution of Sersic indices in core ellipticals is relatively narrow, 
with a characteristic value somewhat higher than those in cusp 
galaxies (as expected if they are indeed re-merger remnants). 

Figure |28] compares the Sersic indices of the outer com- 
ponent in our standard two-component decompositions fitted 
to our gas-rich merger simulations, cusp and core ellipticals. 
Here, the same fitting procedure has been applied over a self- 
consistent dynamic range, regardless of the cusp or core status 
(profile at < 50pc) of the system (for a detailed discussion of 
the dynamic range require ments for reliable Sersic index fits, 
see iKormendv et alllioOSb . As we demonstrated in Paper II, 
for cusp ellipticals and gas-rich merger remnants considered 
in this fashion, there is no dependence of outer Sersic index 
on stellar mass or other galaxy properties. The cusps and gas- 
rich mergers have a relatively narrow Sersic index distribution 
centered around «s ^ 2.6-2.7 with a width of ±0.75. 

Applying our two-component models to the core ellipticals, 
we obtain a similar result. There is no evidence for a signif- 
icant dependence of the fitted outer Sersic index on mass or 
any other property within the core elliptical sample. There is 
an offset between the Sersic index distribution of this sample 
and that of the cusp ellipticals, roughly corresponding to an 
increase in Sersic index of Awj « 1 (i.e. a median '^3.6 
and dispersion ±1.0). 

As described in § ID the change in the outer profile shape 
from cusp to core ellipticals is expected to be a real con- 
sequence of re-mergers. Figure |29] illustrates the difference 
between typical profiles before and after a re-merger. To 
lowest order, the profile is simply puffed out by a uniform 
factor (which has no effect on the Sersic index). However, 
there is significant scatter in final radii at fixed initial radii 
- in other words, more stellar light will be scattered to large 
radii with a re-merger (and, roughly speaking, this fraction 
will increase weakly as the amount of violence, in terms of 
the orbital energy, of the re-mergers is increased; see e.g. 
iBovlan-Kolchin et al.ll2006 ^. This will slightly raise the outer 
Sersic index, and help build up the extended envelopes seen 
around many massive galaxies and BCGs. Exactly how much 
this outer profile is raised depends somewhat on how we de- 
fine our profile fits - in a more robust physical sense, it reflects 
a lognormal scattering of stars by ^ 0.4 dex (Figure [U. 

Since there was no strong dependence on mass in the cusp 
population, we expect no strong dependence in the core popu- 
lation, so long as there is not a dramatic dependence of merger 
history (within this sub-population) on mass. To the extent 
that we see these effects and no substantial mass dependence, 
it implies that most of our core ellipticals have undergone a 
relatively small number of re-mergers, and are therefore rea- 
sonable analogs of our re-merger simulations (where we see 
this change - where we see large envelopes that imply very 
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Fig. 28. — Outer Sersic indices (from our two-component decomposition) in ellipticals and simulated gas-rich merger remnants (point style as in Figure [TSt . 
Filled diamond is typical sightline-to-sightline variation in the simulations, open diamond the source-to-source (or band-to-band) scatter in observed profile fits. 
Gas-rich merger remnants and cusp elHpticals have characteristic «i ~ 2 — 3, without a strong systematic dependence on mass or other properties (dotted line is 
the median for cusp ellipticals and gas-iich merger remnants, = 2.7); core ellipticals are the same, but with slightly higher characteristic h, ~ 3 - 5. 

suits are based on fitting those profiles to either single Sersic 
or core-Sersic type laws (or similarly an outer Sersic law al- 
lowing for some inner "deficit" or core), and are not directly 
comparable to the results fitting the two-component models 
adopted here. For comparison, we therefore consider the re- 
sults obtained with these fitting functions and how they relate 
to the galactic components inferred from our two-component 
decompositions. 

Figures|29]&[30]compare the results of our two-component 
fits with those obtained fitting our simulated re-mergers to a 
"core-Sersic" law of the form 
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Fig. 29. — Illustration of how profile shapes change in a re-merger, as a 
consequence of the scattering of stars. Top: Gas-rich merger remnant. Dot- 
ted blue line shows the outer Sersic index from our typical two-component 
fits (n, = 2.6; concave down in this projection) with the rms deviations of 
the profile from the best-fit two-component model (the rise above this curve 
accurately traces the physical starburst component). Dotted red line shows 
the result of fitting the entire profile to a core-Sersic profile (or a Sersic pro- 
file with some downward break at small radii; ris = 4.1). Both methods yield 
reasonably good fits, but the meaning and values of are different (the n, 
obtained fitting the profile in the latter manner reflects some combination of 
the dissipational and dissipationless components). For gas-rich merger rem- 
nants, the sharper "break" in the profile also yields significant residuals about 
the cored-Sersic fit. Bottom: Same, after a gas-poor re-merger The two- 
component decomposition is just as successful at recovering the physical dis- 
sipational and dissipationless components. However, the scattering of stars to 
larger radii and smoothing of the transition to the starburst component yields 
a situation where other possible fitting functions - in this case a core-Sersic 
fit, with a higher outer index (n, = 5.2; concave up here) and central deficit - 
can be excellent fits in a formal statistical sense. 



high Sersic indices for the outer components - is at the ex- 
treme masses of > a few L» BCGs, where we have empha- 
sized our models should not be extrapolated without caution). 

This relatively modest increase in Sersic index may seem 
surprising, as are the (still relatively low) Sersic indices fitted 
in this manner, given that massive ellipticals are often fitted to 
high Sersic indices > 6. However, we caution that those re- 



(e.g. iGraham et al]|2003h . where rt, is the core break radius 
within which the profile shifts to a power law of slope 7, re 
is the effective radius and n the Sersic index of the outer light 
profile, and a a parameter describing how rapidly the break 
occurs. Whether we free a or set a ^ 00 (the limit of a sharp 
break into the profile core) we obtain similar results (although 
the fit quality x^/i^ is slightly improved with a free a). 

We also consider fitting a single Sersic index to the outer 
profile, i.e. the pro file outside of the cent ral core (following 
the methodology of lKormendv et al.ll2008h . In practice, we fit 
a core-Sersic law to identify the radius where the inner profile 
breaks from the outer Sersic law, and then re-fit the radii larger 
than ^1-2 times this radius to a single Sersic law (we can 
also in more detail fit steadily inwards from larger radii until 
the fit residuals begin to rapidly grow, which yields similar 
results). 

In the case of the direct gas-rich merger remnants, where 
the breaks between the inner and outer components are sharp 
(and where the outer component Sersic indices are relatively 
low), the results are more or less similar The two-component 
fits not only match the physical decompositions in the simu- 
lations, but they also happen to be excellent fits in a formal 
statistical sense. The functional form of the Sersic distribu- 
tion is such that fitting just the outer profiles to a Sersic func- 
tion - in this case giving low n., - yields a shallow inwards 
extrapolation, similar to that in the progenitor (rts ^ 1) disks 
and therefore the dissipationless component; so the extra light 
again appears as a characteristic central excess (although the 
nature of fitting inwards "as far as possible" yields a mild bias 
towards some underestimation of the extra light content). The 
core-Sersic distribution (or equivalently, forcing a single Ser- 
sic function to fit the entire profile, inner and outer), designed 
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Fig. 30. — Top: Surface brightness profiles from tlie smootli extra liglit simulations shown in Figure]!] including the original gas-rich merger remnant (left), 
a major re-merger remnant (center) and another re-merger with different orbital pai'ameters (right). We show the physical decomposition into pre-starburst and 
starburst stars (tracking the starburst component after re-mergers), with the mass fraction in the starburst and Sersic index directly fitted to the pre-starburst 
component. We compare with our two-component fit, with corresponding extra light mass fraction, outer Sersic index, and rms residuals cr^ (in magarcsec"^), 
and show the fit residuals (A^). As before, the two-component decomposition reliably recovers the physical decomposition of the light. Bottom: Same, but 
fitted with a core-Sersic law profile (red) or a single Sersic law fitted only to radii outside of where the central profile begins to deviate from the Sersic law 
(green). In the original gas-rich merger remnant, these methods all show some extra light, and the core-Sersic model performs poorly. After a re-merger, the 
core-Sersic model formally fits the profile well with a higher Sersic index that reflects the combination of the outer profile and the inner extra light profile. If 
we were to simply take /extra as the integrated light relative to that from extrapolation of the Sersic fit to all radii (instead of using our calibrated two-component 
decomposition), we would infer much less "extra" light in the gas-rich remnant and some "deficit" (missing light) in the re-merger remnants - the starburst 
component is still there (with the mass fraction correctly inferred by our two-component fits), but these parameterized fits reflect it only indirectly (with their 
higher Sersic indices) - inferring the "extra" or dissipational component mass from these alternative fitting functions is not physically motivated, even if they are 
better statistical fits to the light profile. 



to describe ellipticals with higher outer and a downward 
break at small radii, performs reasonably well fitting a higher 
n, Sersic distribution to the entire profile (including the extra 
light), but the relatively sharp break between the components 
yields noticeably larger residuals. 

If re-mergers exactly preserved profile shape, the behavior 
would, of course, be identical. In reality however, a re-merger 
scatters stars out to large radii (in addition to the overall "puff- 
ing up" of the profile; Figure [1]), which raises the character- 
istic Sersic index describing the outer profile shape, and 
smooths any sharp breaks in the inner profile. As a conse- 
quence, we obtain rather different results applying alterna- 
tive fitting functions to re-merger remnants. In Figure |29l 
for example, fitting from the outer radii inwards yields a high 
outer Sersic index ~ 5; the nature of the Sersic distribu- 
tion is such that higher-n, outer profiles extrapolate to rela- 
tively large central densities, so the fitted distribution maps 
onto or attaches smoothly to the high central density that is 
the extra light component. Eventually, the inward extrapola- 
tion of the high-n, profiles rises above the light profile (and 
this will be even more pronounced in the presence of black 
hole scouring flattening the nuclear profile slope), so the pro- 
file appears to "breaks downwards" relative to the outer Ser- 
sic fit. Together with the smoothing of sharp features near the 
transition between the dissipational and dissipationless com- 
ponents, this means that one can - after a re-merger - obtain 
very good statistical profile fits with simple cored-Sersic mod- 
els (fjp < 0.04magarcsec~^), with considerably higher outer 
n,. 



Figure |30] demonstrates this in a more rigorous manner: 
we compare the results of these different fitting methods 
with the known physical decompositions in the simulations. 
This emphasizes an important point: the good fits obtained 
with e.g. cored Sersic models (that do not include an ex- 
plicit "extra" component) do not necessarily mean that the 
extra light (as we physically define it: stars formed origi- 
nally in a dissipational starburst) is not still manifest in the 
light profile. This stellar material survives the dry spheroid- 
spheroid re-mergers we simulate, and the rank order of par- 
ticles in radius is sufficiently preserved that it still con- 
tributes a substantial component of the light within ^kpc 
scales in Figure |30] In fact, the high central densities im- 
plied by the large Sersic indices of the fitted cored Sersic 
models (typical n^ ~ 6-8 in observed cored systems; s ee e.g. 
£rugniel & Sim ien|[T997l; [Trujillo et al. 2002; Ferrares e et al.l 
l2006t IKormendv et a l. TOOS) can only be part of the fit if 
the extra light is included - scattering of stars cannot, by 
Louisville's theorem, raise the central density of the elliptical. 
If there were no dissipational component in core ellipticals, 
then their outer portions would resemble a high-n, envelope, 
but their inner portions (< 2 - 3 Re) would resemble a low-n,, 
low density disk (which would mean, among other things, that 
they could not be well-fit over a significant dynamic range by 
any Sersic function). Note that the apparent "deficit" seen in 
these simulations (i.e. shortfall of the inner light profile rela- 
tive to the inwards extrapolation of the outer Sersic fit) does 
not reflect a physical deficit (since the high-«j outer profile 
is much steeper in the central regions than the profiles of the 
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Fig. 31. — Summary of the difference in fitting single gas-rich merger 
remnants (blue) and re-merger remnants (red). We show the distribution for 
both samples of simulations in relative fit quality of a core-Sersic law and our 
two-component decomposition - i.e. the difference in (t^j (in magarcsec"^) 
between the best-fit core-Sersic law and our two-component fit (note the two 
forms have the same number of degrees of freedom). Negative values mean 
the core-Sersic law is formally a better fit. The gas-rich mergers are poorly 
fitted by a core-Sersic law - they almost all prefer the two-component decom- 
position, or little difference (that being in cases where there is only a small 
fraction of extra light). A large fraction of re-merger remnants are better fit 
by a core-Sersic law. 

actual progenitors); this lends some caution to the approach 
of using such comparisons to infer the mass "scoured" by a 
nuclear binary black hole (which we discuss further in §[TUli. 

This general behavior is quite common; Figure |3T] sum- 
marizes these results for our library of simulations - we plot 
the distribution of relative fit quality (for which we adopt the 
proxy of the difference in cr^ between a core-Sersic fit and 
our two-component decomposition), for both single gas-rich 
merger remnants and re-merger remnants. These examples il- 
lustrate the importance of combining models and observations 
in interpreting the parameters obtained from galaxy surface 
brightness profile fits: these different functions are fitting the 
same profiles, with broadly similar quality (and with varying 
relative quality in populations with different merger histories), 
but the Sersic indices and other parameters obtained from the 
fits have different physical interpretations in the context of our 
simulations. Our proposed decompositions, although phys- 
ically motivated, are not required by the observations: in a 
purely statistical sense, the two-component decomposition is 
not only non-unique as a parameterization, it is not even nec- 
essarily the best-fitting choice in all cases. Our intention here 
is not to argue that the proposed decomposition is a better 
fit in a formal sense to observed profiles, but rather that it is 
useful because it tends, when adopted in numerical experi- 
ments, to recover the correct physical decomposition between 
the dense inner (high surface-brightness) component of the 
galaxy originally formed by dissipation of gas into the cen- 
tral regions of the system, and the less dense outer compo- 
nent originally formed by violent relaxation of stellar disks. 
If we were to interpret the extra or excess light to be that light 
in the real profile above the prediction of the best-fit single 
Sersic or core-Sersic model, in the same manner that we do 
relative to the outer components in our two-component fits, 
we would (incorrectly) infer effectively no extra light in most 
of our simulated re-merger remnants. Clearly, calibrating a 
given method against simulations and models is important for 
assigning it a specific physical meaning. 

Interestingly, we find that only our re-merger remnants are 
better fit (in a purely formal sense) by a core-Sersic law, as 
opposed to a two component model, and we find this in gen- 



eral becomes more true as we repeatedly re-merge systems. 
As a consequence, although the parameters of the core-Sersic 
or single Sersic fit may not be directly analogous to those ob- 
tained from our two-component fits, the fact that a galaxy can 
be formally better fit by one of those functions is predicted 
to be a consequence of dry re-mergers which will smooth the 
profile in the manner described above. Indeed, a tight asso- 
ciation is observed between galaxies with nuclear cores and 
galaxies better fit by a single Sersic or core-Sersic law, despite 
differences in the spat i al scales where the two are determined 
IFerrarese etal] 120061: iKormendv et all l2008h . Our models 
provide a natural explanation for this, which not only recon- 
ciles such profile fits with the survival of dissipational/extra 
light components from cusp elliptical progenitors, but in fact 
requires such components to explain the highest observed Ser- 
sic indices from single Sersic or core-Sersic fits. This lends 
further confidence to our modeling of profile shapes, even at 
the level of smooth profiles apparently better fit by functional 
forms different from the convenient two-component decom- 
position we adopt in this paper. 

Given these different choices that can be made when es- 
timating an observational Sersic index, we compare in Fig- 
ure [32] the distribution of Sersic indices for cusp and core 
ellipticals, obtained with different methods. First, we show 
the distribution of Sersic indices measured from our two- 
component decomposition, in both the observed samples and 
our simulations. We compare the distribution extracted from 
our gas-rich merger simulations to the distribution fitted to the 
obs erved cusp el l iptical p rofiles from lKormendy et alJ (1200 8h 
and iLauer et aP (l2007ah . and compare the distribution ex- 
tracted from our re-merger simulations to that fitted to the ob- 
served core ellipticals from the same samples (with the same 
method applied). Because (as demonstrated in Paper II) the 
fitted outer Sersic indices can be quite sensitive to the photo- 
metric depth and systematics between different observations, 
we restrict ourselves to a robust subsample described in Paper 
II: specifically, those ellipticals for which we have > 3 inde- 
pendent sources of photometry that agree to within 20% on 
the fitted value of Wj. We note that we obtain identical results 
(but with larger scatter) if we consider our entire samples. 

We compare these distributions with what we would obtain 
using two other methods to estimate th e outer Sersic index. 
First, we follow iKormendy et alJ (l2008l) and fit a single outer 
Sersic profile inwards from large radii until the central profile 
begins to break upwards or downwards from this extrapola- 
tion, as described above (the comparison in the core ellipticals 
is similar if we fit the entire profile with a core-Sersic distri- 
bution). Second, we attempt to directly fit just the outer radii 
- we first fit a single r'/^ law to estimate the effective radius, 
then consider only the profile at radii larger than ~ 0.5-2.0/?^ 
(within this range, our results are not dramatically sensitive to 
the exact cut). We then fit a single Sersic law just to these 
points. This limits our observational samples to only those 
systems with a large number of observed points at r ^ Rg 
(where such a fit is possible). 

Applied to the cusp ellipticals, these methods yield simi- 
lar distributions, all with median n, ^2.5 and similar scat- 
ter (±0.7-0.9), as expected from our previous comparisons. 
Applied to core ellipticals, the Sersic index distribution is, as 
expected, more sensitive to the method adopted. Fitting di- 
rectly only to large radii yields a roughly similar distribution 
in to our two-component model, with typical '^3-4 
(perhaps somewhat larger). There is more scatter (disper- 
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Fig. 32. — Distribution of outer Sersic indices in ellipticals, using different methods for determining n,, for cusp ellipticals and gas-rich merger remnant 
simulations (top) and core ellipticals and spheroid-spheroid re-merger simulations (bottom). Solid black lines s how the results for o ur entire sample o f simulations 
(each across ~ 100 sightlines). Colored lines show the results for the combined observed samples of Laue r et al] )2007ah and IKormendv et aT (2008), with 
Poisson error bars. Left : Outer Sersic index determined using our standard two-component decomposition. Center: Outer Sersic index determined following 
IKormendv et al.l j20O8l) or Ferrarese et al. (2006) (fitting an "outer" Sersic index inwards, including an extra component if there is some small-scale rise above 
it or allowing a deficit or core-Sersic fit when there is a core or high outer n,). Right: Outer Sersic index determined by re-fitting only the observed profile data 
at r > to a Sersic law. Cuspy ellipticals have a fairly narrow range of n, ~ 2.7 ±0.8, in good agreement with gas-rich merger simulations, and relatively 
insensitive to the parameterization. Core ellipticals also have a narrow range with a higher n ~ 3.5 — 4± 1.0, in good agreement with re-merger sim ulations, 
when parameterized with our two component fit or fitted directly at large radii - however, the issues in the choice of fitting function described in §[T()]make the 
estimated n,, distribution more sensitive to the parameterization. 



sion Ahj ^ 1 - 1.5) - this relates to two effects. First, our 
restriction to radii r > Rg limits the number of data points 
and dynamic range of the fits, increasing the associated un- 
certainties. Second, the fact (as noted above) that the rise 
in outer rts (in our two-component model) in a re-merger is 
primarily driven at large radii from scattered stars, not by 
any increased density at small radii (which still behaves like 
a lower-n, system). Fitting only the large radii is therefore 
not restrained by the curvature and lack of sharply rising den- 
sity at small r, and is more sensitive to the increased scat- 
tering of stars into the envelope at large radii. Adopting the 
method of fitting a single Sersic profile with a downward 
break from Kormendv et al. (2008) or fitting core-Sersic mod- 
els (jOraham et al. 2003; Ferrarese et al. 2006) gives a higher 
typical Hj ~ 4-5 and scatter Arts ^2-3, reflecting the phys- 
ical behavior outlined above. 

In any case, our simulated distributions and the observed 
Sersic index distributions match well when the same meth- 
ods are applied to determine n,. There is a significant (albeit 
not dramatic) and physically meaningful offset in the Sersic 
indices of gas-rich merger remnants and re-merger remnants 
(again, regardless of methodology used to estimate n ,), which 
appears to be reflected in the Sersic index distributions of cusp 
and core ellipticals. This supports the notion that there is some 
difference in the populations, in the sense expected for a rela- 
tively small number of major dry re-mergers. 



10. THE APPEARANCE OF "MISSING" LIGHT: A CAUTION 

As we discussed in § |4] and § |9] core ellipticals are some- 
times referred to as having "missing light" or "mass deficits" 
at small radii (generally below the size scales of the extra light 
of interest to us here), which has been associated with the 
action of black hole "scouring" scattering stars from the nu- 
cleus. We have argued that such scouring should not make 
a difference to the scales of interest for our inferences of the 
dissipational content in these objects. Having said that, how- 
ever, the co-existence of "extra" and "missing" light (i.e. the 
presence of cores in the central profiles of ellipticals domi- 
nated, at their center, by the dense remnants of dissipational 
star formation) can mean that some additional care is needed 
when observers attempt to estimate the amount of "missing" 
light (i.e. how much stellar mass has been kicked out from 
small radii in order to form the central core). The problem is 
that one desires a method to empirically determine where (and 
by how much) the profile has been modified by scouring. Of 
course, any such estimate requires some knowledge of what 
the profile shape was before the scouring took place (roughly, 
the pre-"dry merger" profile of the galaxy). If the galaxy pro- 
file were simply one-component, then provided some physi- 
cal model for the profile or robust constraints on the profile 
shapes of progenitor (non-scoured) galaxies, this would be 
straightforward. But the profile of the progenitors is in fact 
multi-component - so the question arises, can (and if so, by 
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how much does) the amount of "extra light" or dissipation in 
the galaxy bias the estimators of scoured/"missing" light? 

Our comparisons in § |9] should lend some insight to this 
discussion. In Figure |30] for example, we consider the pro- 
files of re-merger remnants fit to both our two-component de- 
composition and to (instead) a core-Sersic distribution or an 
outer single Sersic distribution down to radii including the 
outer portions of the dissipational/extra light component. If 
we were to assume that the true "pre dry-merger" profile of 
the galaxies was given by the inwards extrapolation of what- 
ever fitting law we chose (or, in the case of the core-Sersic law, 
by the inwards extrapolation of the "Sersic" part of the law, 
without the central flattening implied by the "core"), then the 
application of the core-Sersic law to the simulations in Fig- 
ure |30] yields some apparent "missing light" in their nuclei 
(i.e. they fall below the inward extrapolation of the outer Ser- 
sic indices from these particular fitting functions). This hap- 
pens not because of scouring in our simulations (recall, we 
do not resolve the scales where initial cusps would be trans- 
formed into cores), but because the outer Sersic indices of the 
core-Sersic (or limited single Sersic) fits are quite large, and 
large Sersic indices extrapolate (as r ^ 0) to very high central 
densities - even above the true densities of the dissipational 
central components (which themselves drive the fits towards 
high n.5 values, as described in § |9]l. In short, the estimate 
of the "pre-merger" (or "un-cored") profile at small radii is 
biased to higher values than the true "pre-merger" profile. 

Physically, two things contribute to this bias (see § 
the re-merger smooths the difference between the extra light 
(which rises to higher surface densities at r < Re) and dissipa- 
tionless/outer component, leading to some of it (potentially) 
being included in the fit; also, the re-merger scatters more 
stars out to large radii (by ^ 0.4 dex), broadening the overall 
profile and raising the outer Wj. 

This can yield a bias because of the nature of the Sersic 
profile - if we fix the profile normalization at Re, then chang- 
ing «j from, for example, n., = 3 to = 6 yields relatively 
little change at moderate radii (amounting to a factor ~ 2 - 3 
higher surface brightness at r = Q.lRe and r = 5 Re), but as 
r — *■ the profiles diverge strongly. An = 3 profile contains 

0.1 -0.2% of its mass within 0.01/?,., but an «, = 6 pro- 
file contains ~ 1-2%, an order of magnitude difference. If 
the nuclear profile shape within Q.QlRe were flat in both 
cases, one would therefore infer an order of magnitude more 
"missing light" in the higher-n, case. In fact, all the observed 
cases where the "missing light" fraction is estimated to be as 
high as > 1% (~ 10 - 100 Mbh) have high outer Sersic in- 
dices ^ 6 - 10 (seelFerrarese et al. l l2006t ICdte et al.ll2006t 
iKormendv et al.l l2008lK The physical concern is that these 
Sersic indices are much higher than those of typical "cusp" 
ellipticals (with slopes corresponding to Wj ~ 2 - 4 in their 
central regions) - in other words, the assumed "pre-scouring" 
profile in these cases is rising much more steeply as r ^ 
than is actually observed in the presumed progenitors of core 
ellipticals. 

This is not to say that the presence of central light deficits 
is meaningless - simply that, with this methodology, the in- 
ferred amount of missing light at small scales (^ Rbh) can be 
biased (generally towards higher values) by features on large 
scales Re). A more robust null hypothesis (i.e. estimator 
of what the nuclear profile would be in the absence of scour- 
ing) may be obtained if we assume the "pre-merger" nuclear 
profile rises about as steeply as observed in typical cusp ellip- 
ticals. 



For example, if we examine the profile of just the inner 
dominant component (the "extra light") in cusp ellipticals, we 
expect that the median non-scoured distribution would rise 
about as steeply as an Hs ^3.6 profile (or a Nuker profile with 
inner logarithmic slope 7 0.7-0.8) once well interior to the 
effective radius of the extra light. In core ellipticals, the ob- 
served nuclear profile at r ^ Rexti^ is flatter, more like median 
n.s ^1.6 (or Nuker slope 7 ~ 0.2) - this, of course, is why 
they are cores. Taking the typical nuclear shape in cusp ellip- 
ticals as our null hypothesis, we obtain characteristic missing 
light fractions ~ 0.5-3 x 10~-'Mc,ai, comparable to the black 
hole mass and expectation s from scouring models. Alterna- 
tively, Lauer et aLl (l2007bl) estimate the scoured mass as that 
within (defined as the radius where the logarithmic slope 
7 = 1/2, which is reasonably representative of the radius in 
which scouring has significantly flattened the profile), giving 
a median value ^ lAM^u- 

This is important because the high values sometimes in- 
ferred for the scoured mass fraction or "missing light" (> 
0.01 Mgai) would actually present a challenge to the black hole 
scouring models - those models predict a scoured core mass 
~ Mbh, but in some cases these empirical estimators yield 
scoured masses 10- 100 Mbh- Our simulations demonstrate 
that these extreme observational estimates can probably be 
explained as artifacts of the fitting functions chosen, in cases 
where the outer «, is large. Qualitatively, this is not expected 
to alter the classification of cusp and core ellipticals; but if 
a robust quantitative estimate of the scoured or core mass is 
desired, then care is needed in empirically evaluating these 
cases. We find that if we adopt estimators less sensitive to the 
profile at large radii, these high values come down, more in 
line with the expectations from scouring models. 

11. DISCUSSION 

In Paper I and Paper II, we demonstrated that "extra light," 
in the sense of stellar populations formed from dissipation in 
e.g. a gas-rich merger induced starburst, is ubiquitous in cusp 
ellipticals and gas-rich merger remnants, and studied how its 
properties are related to and (in some cases) drive those of 
the galaxy. Here, we show that this dissipational component 
is expected to survive subsequent re-mergers of those ellipti- 
cals, even major gas-poor spheroid-spheroid mergers, in the 
sense that it will continue to contribute substantially to the 
central light profile and can be empirically recovered. We 
apply this to a large observed sample of ellipticals with cen- 
tral cores (i.e. flattening of their light profiles within the cen- 
tral ^ 30-50pc), and show that they are consistent with a 
surviving dissipational "relic extra light" component which 
our adopted empirical fitting machinery and comparison with 
simulations allow us to recover and interpret. 

11.1. Comparing Simulations and Observations: Empirical 
Decomposition of Light Profiles 

In Paper II, we argued that stars in cuspy ellipticals/gas- 
rich merger remnants should be separated into at least two 
distinct populations. First, stars which are formed in the disks 
(or otherwise in extended distributions in progenitor galaxies) 
before the final merger and coalescence of the progenitors. 
The final merger scatters these stellar orbits and they undergo 
violent relaxation. They dominate the light, even in highly 
gas-rich merger remnants, outside of 0.5 - 1 kpc, and form 
a Sersic-law profile owing to their partial violent relaxation. 
Second, the starburst or dissipational population, formed in 
the central gas concentration in the final merger This compo- 
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nent is compact, and dominates the light inside a small radius 
< 0.5- Ikpc. These stars do not undergo significant violent 
relaxation, but form in a nearly fixed background potential set 
by the dissipationless component of the galaxy. We devel- 
oped, tested, and studied in detail in Paper I a two-component 
fit decomposition that, in simulations, could reliably extract 
the properties of these two physically distinct components. 

Here, we demonstrate that, to lowest order, the two pri- 
mary components structuring the surface brightness profile 
are expected to survive in dissipationless spheroid-spheroid 
re-mergers, in the sense that the physical starburst component 
still forms a more compact distribution that dominates or con- 
tributes significantly to the profile at small radii. Re-mergers 
will generally expand or puff up the system by a factor ^ 2 
in size, and smooth out the profile by mixing stellar popula- 
tions and scattering stars by a typical factor ^ 0.4dex in radius 
(see Figure [T]) - this can make fitting the profile more sensi- 
tive to the prescription adopted, and may smooth an obvious 
break in the profile around the transition from dissipational 
to dissipationless components (Figures |29ll30ll, but does not 
fundamentally remove the extra light in a physical sense. 

We apply our same fitting procedure to simulated re-merger 
remnants, and find that, despite the re-merger, it is still able to 
statistically recover the physically meaningful components of 
the original, spheroid-forming gas-rich merger (i.e. the origi- 
nal dissipational or starburst component and the original dis- 
sipationless or pre-starburst component; see Figures |2]|4l). In 
other words, even after a re-merger, the surface brightness 
profiles of ellipticals retain information about the gas frac- 
tions and starburst mass fractions of their original gas-rich 
mergers. Our parametric fitting form is simple. We consider 
the sum of two Sersic laws: an inner extra light component, 
for which a fixed = 1 works best in a mean sense when the 
data are not especially constraining in the central regions, but 
for which can be freed if the data are of sufficient quality; 
and an outer component with a free Sersic index. We explic- 
itly demonstrate that this approach is successful even when 
(as is common in re-merger remnants) the profile is smoothed 
by re-mergers and obvious breaks might not be present in the 
profile (FiguresEllEB. 

We apply this to a large sample of "core" ellipticals, which 
(it is often argued) have been modified by re-mergers, and 
find that it is reliable and implies significant "extra light" (by 
which we mean the remains of the original gas-rich compo- 
nent) in almost all cases (Figures [SlfTTb. We also match each 
of the observed profiles to our library of simulations (of both 
gas-rich mergers and re-mergers) - i.e. directly find the sim- 
ulation mass profile which most closely resembles that ob- 
served, and consider non-parametric estimators of the mass 
in a central starburst component. We find that in all cases we 
have simulations which provide good matches to the observed 
systems, comparable to the typical point-to-point variance in- 
herent in the simulation surface brightness profiles (cr^ ;$ 0. 1). 
The physical starburst components in these best-fitting simu- 
lations are closely related to those that we fit directly to the 
observed profiles, lending further support to our attempt to 
physically decompose the profiles (Figure [TsTl. Where avail- 
able, stellar population models including extended star for- 

There is also a third component present in simulations and impor- 
tant for observed kinematics but not prominent in light profile fitting: gas 
moved to large radii temporarily either by feedback or tidal effects, which 
settles into the relaxed remnant and re-forms small rotationally supported 
components (embedded disks, kinematically decoupled cores, etc.; e.g. 
IHemquist & Bamesll99ltlHopkins et aU2009l) . 



mation and a subsequent burst independently support our in- 
ferred starburst mass fractions (Figure [TSll. Likewise, metal- 
licity, age, and abundance gradients, where available, support 
our decompositions, typically demonstrating a smooth tran- 
sition to a younger, more metal rich population at the radii 
where the dissipational component begins to dominate the 
profile, as predicted in Paper II (Figure [T9]l. A complete list 
of fit parameters and compiled galaxy properties is included 
in Tableffl 

11. 2. Predictions and Observations 

Fundamentally, we argue that all ellipticals - including 
those with central cores - are "extra light" ellipticals (dis- 
sipational systems), insofar as "extra light" refers to a com- 
ponent originally formed dissipationally on top of a back- 
ground of dissipationlessly scattered envelope (a physical 
two-component nature). Ellipticals with cores show just as 
much dissipation, at a given mass, as cuspy or power-law el- 
lipticals. Over the core population, the mass contributed by 
a central starburst or dissipational population can vary sub- 
stantially, but in a physical sense, it always exists in addition 
to the pre-starburst or violently relaxed stellar populations. 
In terms of their dissipational or starburst population prop- 
erties, we demonstrate that core ellipticals appear to be, for 
the most part, a continuous extension of the cusp population. 
The correlations obeyed by the inferred dissipational compo- 
nent itself are continuous and agree where there is overlap in 
galaxy properties between the cusp and core populations, as 
expected in physical models where core ellipticals are formed 
by re-mergers of cusp elliptical progenitors. 

Specifically, we find: 

(1) The mass fraction in the dissipational or starburst com- 
ponent of both cusp and core ellipticals is a strongly decreas- 
ing function of mass (Figures I20II2TI 1. In detail, the mean star- 
burst mass fractions can be approximated as Equation (O: 

with a factor ^ 2 scatter at each mass. The trend is similar for 
both cusp and core ellipticals, and the gas fractions needed 
span a range bracketed by the typical observed gas fractions 
of spkal galaxies at the same mass, at z = (bracketing 
the low end of the required gas fractions) and z ~ 2 - 3 
(bracketing the high end). Core ellipticals have not only 
preserved their dissipational content, but its fraction reflects 
that of their ultimate progenitor disks. 

(2) At each mass, the degree of dissipation strongly affects 
the sizes of the remnants. In both observations and simu- 
lations we demonstrate a tight correlation between effective 
(half-light) radius at a given stellar mass and the inferred 
dissipational/extra light fraction (Figure |24] |. This owes to 
the compact nature of the central dissipational component - 
increasing the mass fraction in this component means that the 
half-light radius must be smaller The correlations obeyed be- 
tween dissipational content and size at fixed mass are similar 
for cusp and core ellipticals. 

Despite the fact that re-mergers puff up remnants by a 
factor ^ 2, their masses also double, so this only moves them 
by about ^ 0.1 dex off the mean size-stellar mass relation 
obeyed by gas-rich merger remnants (which has an intrinsic 
scatter ~ 0.3 dex). On the other hand, we demonstrate that 
changes in the original starburst or extra light fraction can 
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alter the remnant size by nearly an order of magnitude at 
fixed mass. Thus, even if core ellipticals have undergone a 
moderate number of re-mergers, we expect that the degree of 
dissipation in their original formation should still be the most 
important factor setting their sizes today, and we demonstrate 
this in the observations. 

(3) Re-mergers roughly preserve the size-mass relation of 
the dissipational/starburst component itself, approximately 
set by the radius at which the starburst component becomes 
self-gravitating (GMexna/^extia ~ GM,//?(.); and we do find 
that the core ellipticals and cusp ellipticals obey nearly 
the same correlation (Figure |22] |. If re-mergers preferen- 
tially puff out low binding energy material (i.e. material at 
large initial radii), then we might expect that the quantity 
(GMextra/^extra ~ should increase slightly in a re- 
merger (by a maximum factor ~ 2 in a 1:1 merger, if the inner 
component is not expanded at all, and the outer component 
expands by a factor of 2). We see tentative evidence for such 
an offset, but it is small (a factor ^ 1.4- 1.5; Figure |23]|. 

(4) In Paper II we demonstrated that the extra light compo- 
nent gives rise to stellar population gradients in the remnant. 
We show here that these gradients are only weakly affected 
by re-mergers, in agreement with observations which find 
that cusp and core ellipticals do not show significant differ- 
ences (at fixed galaxy properties) in their gradient strengths 
(Figure |6]|7]l. To the extent that both radius and gradient 
strength scale with the contribution of the original dissipa- 
tional component at fixed mass, we find that our prediction 
from Paper II should also hold for core ellipticals (namely 
that at fixed mass, smaller ellipticals should, on average, 
exhibit stronger metallicity gradients), and e arly observa- 
tional evidence appears to support t his (e.g. iMehlert et alJ 
12003'; Sanc hez-Blazquez et al.1 l2007t iReda et all l2007 l see 
also Figure [T9] l. Ther e are other Unes of evidence for the 
survival of extra light: iLauer et aH (l2005h find that nuclear 
color gradients, ellipticities, and isophotal twists in the 
centers of cusp and core galaxies are continuous and trace 
similar distributions as a function of mass and luminosity. 
Analogous correlations should hold with integrated stellar 
populations, which we discuss in Paper II, provided that these 
are correlated with the gas fractions of the pre-merger disks. 

(5) Given the appropriate dissipational mass fraction 
inferred from fits to the observed surface brightness profiles, 
simulated re-merger remnants reproduce the global kine- 
matics and isophotal shape distributions of core ellipticals 
(specifically rotation V/cr, ellipticity, and boxyness 04/0), 
and their mutual correlations (Figures I26ll27] i. These distri- 
butions are only reproduced by simulations with matched 
dissipational mass fractions: systems with too little dissi- 
pation (lacking a high central density to make the potential 
more round and disrupt box orbits) are too elliptical, and 
systems with too much dissipation remain too disky and 
rapidly rotating even after several re-mergers (Figures |25] & 
|27] |. More detailed comparison of simulations with full two- 
dimensional velocity fields, along with higher-order measures 
of the velocity field at a given point, are needed to robustly 
discern whether scenarios as simple as a single re-merger 
are viable formation mechanisms for slowly-rotating, boxy 
ellipticals, but similar dissipational content appears to be a 
basic prerequisite. 



(6) We predict and find that the outer Sersic indices of el- 
lipticals with cores (in our two-component decompositions) 
trace a roughly constant distribution with median ^ 3-4 
and scatter An, ^ 1, without a strong dependence on galaxy 
mass, effective radius, or other properties (Figure |28] ). We 
found similar results in Paper II for the cusp elliptical pop- 
ulation, but with a lower median ^ 2-3. The difference 
arises naturally in our re-merger simulations, as the scatter- 
ing of stars in subsequent mergers will tend to broaden the 
outer regions of the light profile, leaving a less steep falloff 
and giving rise to a higher Sersic index. That there is an off- 
set in the outer Sersic index distribution of these populations 
suggests that indeed some sizable fraction of core ellipticals 
have experienced a re-merger, but that the offset is relatively 
small suggests that there have not been a large number of such 
re-mergers. 

We emphasize that these Sersic indices are not directly 
comparable to those in previous studies, which fit different 
functional forms to the light profile or fit e.g. a Sersic pro- 
file to the entire light profile (including the extra light). If 
we fit our simulations to a single Sersic index or core-Sersic 
profile where it is the formal best fit, we recover a stronger 
dependence of Sersic index on galaxy ma ss, luminosity, or 



effective radius, similar to earlier claims dCaon et al 
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'Prugniel & Simien 1997; 'Graham"200l iTruiillo et al 
Ferrarese et al. 2006): our models are consistent with the 
results derived using these methods, however this does not 
strictly trace the physical dissipationless component of the 
galaxy. In fact, we show in §l9](Figures l29ll3TI ) that re-merger 
remnants are often formally well-fit by core-Sersic laws (Ser- 
sic profiles with a central deficit), in agreement with obser- 
vations of core ellipticals. These cases obtain high Sersic in- 
dices, driven in part by the presence of a dissipational com- 
ponent at small radii, raising the central surface brightness 
(Figure [29]). These functional forms are a precise formal pa- 
rameterization of the light profile, but the best-fit parameters 
determined in this manner reflect the combination of the dis- 
sipational and dissipationless outer components, and do not 
trivially translate to physical descriptions of the components 
of the galaxy. 

Indeed, as demonstrated in Paper II, many of the results 
in studies using these fitting approaches are actually driven 
by a dependence of extra light on galaxy properties. The 
meaningful dependence of outer Sersic index as we quantify 
it (tracing the dissipationless component) on mass owes to 
the dependence of typical merger history on mass (dwarf 
spheroidals and pseudobulges will have lower typical 
than gas-rich merger remnants, which themselves have lower 
typical Hs than re-merger remnants, and these populations, on 
average form an increasing sequence in mass). For a broadly 
similar merger history (in terms of e.g. number of major 
mergers), the profile of the true dissipationless component is 
expected to be self-similar, as we recover (for both the cusp 
and core populations), since it is determined purely by gravity. 

These and other correlations argue that core ellipticals 
should be properly thought of as containing "extra" or dissipa- 
tional components in the same physical (although not neces- 
sarily observational) manner as cusp ellipticals; their central, 
high density regions formed via dissipational processes. We 
therefore refer to Paper II for a number of other proposed ob- 
servational tests of the models herein, which should hold for 
core ellipticals as well. This is a critical test of the merger 
hypothesis and supports the notion that even massive, slowly 
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rotating, boxy ellipticals with cores were originally formed in 
gas-rich major mergers, albeit potentially modified by subse- 
quent re-mergers. 

Points (1) and (2) emphasize that, even if a system has ex- 
perienced re-mergers, and even in observed massive core el- 
lipticals, dissipation is the key driver of galaxies along the 
fundamental plane (in terms of galaxy stellar mass). The de- 
gree of dissipation in the original gas-rich merger is the most 
important factor determining the size of the object and the ra- 
tio of baryonic to dark matter within the effective radius of 
the stellar light, and the correlations obeyed by the observed 
systems agree well with our simulations. 

We show that this degree of dissipation is a systematic func- 
tion of mass, without a significant offset between cusp and 
core elliptical populations (as expected if it is set by the gas 
fractions of the progenitor disks, and there have not been a 
large number of re-mergers of individual objects). In other 
words, we are able to demonstrate that the amount of dissipa- 
tion expected based on known disk gas fractions as a function 
of mass is precisely that needed to explain the extra light in 
the surface brightness profiles and to reconcile the densities 
and radii of disks and e llipticals as a fu nction of mass. We in- 
vestigate this further in lHopkins et al.l (l2008b ). and show how 
it gives rise to the fundamental plane scalings and "tilt." Once 
on the fundamental plane, our simulations and other exper- 
iments (Boylan-Kolchin et al. 2005; Robertson et al. 2006a) 
find that re-mergers typically move more or less parallel to 
the fundamental plane, emphasizing that the tilt must arise in 
the initial, gas-rich spheroid-forming merger 

Although some of these effects will become increasingly 
scattered or smeared out by a large number of re-mergers, 
multiple lines of evidence above suggest that the typical core 
elliptical has undergone relatively limited dry re-merging. If 
re-mergers are indeed the mechanism by which cores are 
formed, our comparisons suggest a small number ^ 1 major 
re-mergers per object since its formation in a gas-rich merger 

This is in line with other observational indications from 
e.g. the mass function evo lution of ellipticals an d red galax- 
ies (e g.lBundv etal.ll2005 ; Borch et al. 2006; Pa nnella et al] 
|2006'; Franceschini et al. 2006; Fontana et al. 2006) and direct 
observational estimates of the dry merger rate (Lin et al. 20()3, 
l2008t Ivan Dokkumll2005l: IBell et alj|2006l). It is also gener- 
ally expecte d in cos mological rnodels (iDe Lucia et al.ll2006t 
IZheng et aL 2007; H opkins et al.l2008^ . where only the most 
massive BCGs withM, > IO'^Mq or so are expected to have 
a large number of major re-mergers. Even in these cases, 
it is not clear whether such mergers actually proceed effi- 
ciently, or whether the secondary is tidally destroyed before 
the merger a nd added as part of an extended envelope or inter- 
cluster ligh t (Gallagher & Ostriker 1972^ iMonaco et aD2006l; 
IConroy et a l. 2007; Purcell et al. 2007). At these masses, the 
growth of the system may also be dominated by a large num- 
ber of minor mergers, rather than a few major mergers (major 
merger s being the dominant "dry" growth mo de at < a few 
L,; see lMaller et ani2006l; iMasjedi et al.ll2008 l). 

Our constraints on this extreme of the population are, unfor- 
tunately, much weaker. Moreover, we do not intend our mod- 
eling to be extrapolated to this regime. Such systems formed, 
for the most part, at early times from potentially very different 
progenitors than those we model, with more complex subse- 
que nt merger hi stories our simulations do not capture (see, 
e.g. i Li et al.ir200 7). More detailed analysis of individual ob- 
jects, and more detailed models incorporating a fully cosmo- 
logical context are called for to study the formation history 



of galaxies in this regime. However, such systems constitute 
only a small fraction (< 10-20%) of the mass density in even 
the core elliptical population (and just ~ 3-5% of the total 
mass density in ellipticals), so the vast majority of the stellar 
populations of spheroids should be reasonably represented by 
our modeling. 

U.3. Profile Fitting and Tests of Core Creation Models 

The existence of nuclear cores in these ellipticals is com- 
monly attributed to "scouring" by a binary black hole in gas- 
poor re-mergers. This will eject stars in close encounters 
with the binary, flattening the initial steeply-rising power-law 
cusps into the observed shallow cores. To the extent that the 
stars in this initial cusp are ejected to larger orbits, leaving 
a core, these remnants can be thought of as having "miss- 
ing" light in their centers. We emphasize that this is com- 
pletely consistent with there also being "extra light," as we 
define it, in these galaxies. The "extra light" we refer to is the 
remnant of stars formed in dissipational central star forma- 
tion events, and dominates the profile within rather large radii 
~kpc, blending smoothly onto the outer dissipationless pro- 
file. Scouring will flatten the nuclear peak or cusp of this dis- 
sipational remnant, but will not be likely to remove the ~ 10% 
of the stellar mass that constitutes the dissipational starburst 
relic. Despite the terminology, "extra" and "missing" light 
are not mutually exclusive. In a strict physical sense, all el- 
lipticals are dissipational/"extra light" ellipticals - they have 
all experienced some dissipational star formation. Of these 
dissipational ellipticals, there are "un-cored" or "un-scoured" 
("power-law," "cusp," or "cuspy core") ellipticals (presum- 
ably objects whose last major merger was gas-rich, allow- 
ing the black holes to merge quickly, and forming the cen- 
tral cusp via new star formation) and "cored" or "scoured" 
ellipticals (presumably the remnants of subsequent spheroid- 
spheroid "dry" re-mergers). The "missing light" is a core - 
a flattened nuclear profile as opposed to a continued steeply 
rising cusp - within the very center of the dense dissipational 
component. 

Our methodology is robust to re-mergers and scouring: our 
two-component profiles still recover the dissipational compo- 
nent in simulations allowing for any reasonable model of nu- 
clear scouring. However, we demonstrate that the same pro- 
file can be well-fitted by different assumed functional forms, 
and the parameters obtained from these fits - particularly in 
regards to the "extra" and "missing" light - are not trivially 
related and require careful interpretation. In particular, obser- 
vations find that core ellipticals can be reasonably well-fit in 
a formal statistical sense by cored Sersic laws (an outer Ser- 
sic profile with an inner "flattening" or deficit), generally with 
high Sersic indices n, > 5, and we obtain consistent results fit- 
ting our re-merger simulations with these methods. There is 
no explicit "extra" component in these fits: rather, in our sim- 
ulations, we find that the high n, values - which extrapolate to 
high central densities before the flattening or deficit - implic- 
itly reflect the extra light (that being the central, dense compo- 
nent that enables the high-n, fit; otherwise the system would 
have a high outer n, but a low central density, not resembling 
any Sersic profile). Because the re-merger has smeared out the 
light profile, removing characteristic features (or kinematic 
subsystems) that might have (in the original gas-rich merger 
remnant) been a more obvious indication of the transition to 
radii where stars formed dissipationally dominate the profile, 
the physical breakdown of these systems is not obvious from 
a strictly empirical standpoint. By calibrating different fitting 



Dissipational Signatures in Core Ellipticals 



39 



methods with our simulations, we have attempting to provide 
an interpretive context and physical motivation for specific in- 
terpretations of these fits. 

In a related manner, this raises some cautions regarding the 
purely empirical methods sometimes used to estimate how a 
profile has been modified by core scouring - i.e. how much 
mass is "missing" from the nuclear region (relative to the 
steep nuclear rise of the progenitor cusps). This is a subtle 
issue. There may of course be no light actually missing after 
scouring - rather, stars are scattered from small radii to large, 
flattening the central profile (what is really desired is an esti- 
mate of how much stellar mass has been scattered - i.e. how 
much stellar mass must be moved from the nuclear region to 
larger radii to explain the difference in profiles). So any such 
estimate is sensitive to the model for the nuclear profiles of 
the progenitors - equivalently, what the profile would be in 
the absence of scouring. 

We explicitly consider one such (commonly adopted) 
model, where the implicitly assumed progenitor profile is 
based on the inwards extrapolation of an outer Sersic law (as 
in e.g. core-Sersic fits). In these fits, the effects described 
above (effects on the profile at larger radii than the core it- 
self) can lead to a very large outer Sersic index in the formal 
best fit; in the high-«j regime (n^ > 6), the nature of the Ser- 
sic profile is such that these profiles then rise steeply at small 
radii (more steeply even than the power-law nuclear profiles 
of cusp ellipticals), and a small change in the outer profile 
(raising n,) can substantially raise the extrapolated (and as- 
sumed "pre-merger") inner profile. This can bias the esti- 
mate of the "missing light" fraction towards very high values 
(~ 1 -5%) and literally interpreted can make it appear as if 
the "core" extends to radii ^kpc. In fact, we find that adopt- 
ing this methodology in such cases, we would infer "missing 
light" at this level even in simulations where there is no phys- 
ical scouring (i.e. there is no actual "missing light"). More- 
over, if real, these values would imply scoured masses ^ 
10-50Mbh, in troubling disagreement with models of scour- 
ing which predict sco ured mass deficits ~ 0. 1-1 Mbh Pgr ma- 
jor re-merger (see e .g. lMilosavlievic et al.l2002l;lMerrittl2006l; 
ISesana et al.l l2007l) . If stars from centrophilic orbits in tri- 
axial potentials allow the binary to coalesce rapidly (provid- 
ing stars from near ^ Rg, where the loss of even IOMbh 
makes no difference to the profile, to ha rden the binary), as 
suggested in some idealized calculations (iBerczik et al.ll2006t 
iHoIlev-Bock elmann & Sigurdsson'20061) and merger remnant 
simulations (Hoffman et al. 2009), the difficulty for the mod- 
els in explaining such large "missing masses" grows. 

Given the steep dependence of the implicitly assumed cen- 
tral progenitor profile on Sersic index in this regime, and the 
appearance of "missing light" even in simulations without 
scouring given a literal interpretation of this fitting methodol- 
ogy, it is likely that the "problems" for the scouring models in 
these extreme cases reflect more the observational uncertainty 
regarding the appropriate "un-scoured" profile, rather than 
fundamental uncertainties in the physics of scouring. Future 
work, as observations and models improve, should attempt to 
carefully model the progenitor profiles and develop estima- 
tors that are less sensitive to the profile shape at radii much 
larger than the core itself. An important test for any such es- 
timators should be that re-merged progenitors with initially 
cuspy profiles (down to the desired resolution limits), without 
scouring included, should typically yield little or no "miss- 
ing light." We briefly consider a couple of possibilities for 
such estimators of the scoured mass (less sensitive to the large 



scale radii), and obtain (re-analyzing the observations) typical 
scoured mass estimates of ~ 0.5 - 3 Mbh, in better agreement 
with scouring models. 

Further study, in particular observation of the nuclear re- 
gions ^ 7?BH is needed to test models of scouring and test the 
accuracy of different estimators for the scoured stellar mass. 
For example, scouring is expected to preferentially eliminate 
stars on radial orbits and leave a bias for ta ngential orbits 
within the radius aff ected (e.g. lOuinlan & H ernquist 1997). 
iGebhardt et al.l (|2003) see tentative evidence for this in a lim- 
ited sample of ellipticals; the major-axis radii within which 
the effect appears are generally ^ 0.5-37?bh, as expected in 
scouring theories and our revised inferences from core pro- 
files. A preliminary comparison supports our conclusions 
(and caveats) here, but the number of relevant ellipticals is 
small. If similar observations can be obtained for a sample of 
ellipticals with alternative estimates of the missing mass frac- 
tion, it can be determined whether some estimators are biased 
towards putting too much or too little of the profile into the 
scoured component. 

11.4. Summary 

We have developed a paradigm to understand the structure 
of both cusp and core ellipticals, in which there are funda- 
mentally two stellar components: a dissipational central star- 
burst component and a more extended violently relaxed com- 
ponent. We have shown that the separation between these 
components can be inferred with observations of sufficient 
quality, and used to understand the formation history of ellip- 
ticals as a function of a wide range of properties. This allows 
us to demonstrate that dissipation is critical to understanding 
the properties of ellipticals, including (but not limited to) the 
structure of their surface brightness profiles, their sizes, ellip- 
ticities, isophotal shapes and rotation, age, color, and metal- 
licity gradients (and their evolution), and the gas content and 
properties of their progenitors. 

In particular, we argue here that this remains true for ellip- 
ticals with cores and re-merger remnants - in other words, all 
ellipticals, including core ellipticals, are fundamentally "ex- 
tra light" ellipticals. The core ellipticals form a continuous 
family with the cusp ellipticals in terms of their dissipational 
content, as expected if these cusp ellipticals are their progen- 
itors and if the number of re-mergers (if important) has not 
been large ('^ 1 - 3 major re-mergers) for the typical core el- 
liptical. We demonstrate that, despite the possibility that these 
systems have expanded via re-mergers, the degree of dissipa- 
tion in the original gas-rich merger (the memory of which is 
retained in the surface brightness profile) remains the most 
important factor determining the size, gradient strength, and 
other properties of the remnant. 

We have studied the properties and identified robust trends 
of dissipational stellar remnants in the nuclei of elliptical 
galaxies with cores and remnants of gas-poor, spheroid- 
spheroid re-mergers, across a large set of simulations, in 
which we vary e.g. the galaxy masses, initial gas fractions, 
concentrations, halo masses, presence or absence of bulges, 
presence or absence of black holes, feedback parameters from 
supernovae and stellar winds, orbital parameters and disk in- 
clinations, and mass ratios of the merging galaxies. This 
range of parameters allows us to identify the most important 
physics. 

As we found in Paper I and Paper II, the most important 
factor determining the structure of the remnant (insofar as the 
properties we consider are concerned) is how much mass is in 
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the original (in the original or last spheroid forming, gas-rich 
merger) dissipationless (violently relaxed) component, versus 
the mass fraction in the dissipational (starburst) component. 
Orbital parameters and initial galaxy structure can, in princi- 
ple, affect the remnant surface brightness profile significantly, 
but only indirectly, insofar as these help to set the amount of 
gas which will be available at the time of the final coalescence 
of the galaxy nuclei (i.e. how much mass ends up in the star- 
burst component, as opposed to being violently relaxed in this 
final merger). 

Re-mergers will expand the original gas-rich merger rem- 
nants, and smooth the profile (scattering stars substantially 
about their mean final radii), but, as it is well established 
that dissipationless mergers conserve (in a mean sense) parti- 
cle rank order in binding energy and therefore radius (Barnesj 
p88), they will preserve these components. In other words, 
the product of a re-merger is, to lowest order, the sum of the 
two progenitor spheroid dissipationless components (consti- 
tuting the nearly self-similar outer violently relaxed compo- 
nent) and their inner dissipational/starburst components (con- 
stituting the central dissipational/starburst remnant compo- 
nent of the final re-merger remnant, despite the fact that there 
might be little or no new dissipation in the re-merger). 

We have demonstrated that this makes predictions for how 
fundamental plane scalings arise, which we study further in 
[Hopkins et al. (2008b). We make a wide range of new pre- 
dictions for the distributions of these properties and how they 
scale with the degree of dissipation, and how they should scale 
with each other and various other observational proxies for 
this degree of dissipation (which we define herein). We have 
predicted and shown (given these proxies) that dissipation is 
indeed more important (contributing a larger mass fraction) in 
low-mass ellipticals, in line with expectations based on how 
gas fractions are known to scale with disk mass. Testing all of 
these with better observations should be possible in the near 
future, with well-defined samples of ellipticals and continued 
improvements in mapping e.g. the surface brightness profiles, 
stellar populations and their gradients, and structural proper- 
ties of ellipticals over a wide dynamic range. 

Given our decompositions, we observationally confirm the 
long-standing prediction of the merger hypothesis, that suf- 
ficient dissipation should have occurred in the inner regions 
of ellipticals to explain the discrepancy between their central 
densities and those of their progenitor spirals, a confirmation 
that fits well in line with what is now well-established in gas- 
rich merger simulations and is also directly seen in progress 
in ongoing/recent mergers, which have (through clear re- 
cent central star formation) raised their phase space de nsities 
to be comparable to ellipticals jKo rmendy & Sander ^ [ 19921: 
Dovon et a lj Il994t iGenzel et al/ 2001: Ro thberg & JosephI 
2006; Tacc oni et al1 [2002: Dasyra et al. 200?!). We show that 
this is true for the core elliptical population just as we demon- 
strated it for merger remnants in Paper I and cusp ellipticals 
in Paper II. In other words, the same mechanism can explain 
the central densities in both cusp and core ellipticals - while 
core ellipticals may be modified by subsequent re-mergers, 
no alternative formation mechanism for them (in the sense of 
formation of their spheroid progenitors) is required. 

We demonstrate some important caveats in observational 
studies of core ellipticals. These ellipticals are multi- 
component (dissipational plus dissipationless, in the physical 
sense above) galaxies, but they have shallow nuclear profiles 
(central cores) instead of steeply rising central cusps, and can 
therefore be thought of as "missing light" galaxies in the sense 



that some process (e.g. scattering of stars from a black hole bi- 
nary) has flattened the nuclear profile slope to create the core. 
This is fundamentally different from the processes involved 
in the evolution of the "extra light" (by which we mean the 
remnant of the starburst/dissipational star formation) and they 
should not be confused - whatever process forms cores likely 
acts on small scales and involves a relatively small fraction of 
the galaxy mass, much less than the '^kpc scales and ~ 10% 
of the galaxy mass that are characteristic of the dissipational 
component. Owing to the smoothing of the outer profile in 
a re-merger, re-merger remnants can often be better fit (in a 
purely formal sense) by e.g. single Sersic profiles with central 
deficits or core-Sersic laws: where observed, we argue that 
this is an indication of such re-mergers, but we stress that it 
does not mean the galaxies are not, in fact, two-component 
objects in a physical sense (simply that the fitted function re- 
flects some combination of the two components). For this 
reason, the two-component profiles adopted here lend them- 
selves to more direct physical interpretation and comparison 
with galaxy formation models. Furthermore, we demonstrate 
that care is needed when using such core-Sersic fits as an es- 
timator of the "scoured" nuclear mass - when the outer Sersic 
index is large, the inferred "missing mass" can be quite sensi- 
tive to the details of the profile at large radii. This may explain 
estimates of scoured masses which are much larger than those 
predicted by models of core creation. 

There are, of course, other changes to galaxy properties 
in re-mergers. As discussed in Paper I and Paper II, gas 
which survives the original spheroid-forming gas-rich merger 
will quiescently settle into the galaxy and form kinematic 
sub-components (in particular, embedded disks and kinemat- 
ically decoupled nuclear components). We demonstrated that 
these do not contribute significantly to the surface brightness 
profile, an d therefore they are not evident in our analysis. 
However, Cox et al.' (2006b) and other numerical studies of 
gas-rich mergers (Naab et al. 2006a: Burkert et al. 2008) have 
demonstrated that these subsystems can contribute strongly to 
the isophotal shapes (in particular driving the diskyness of the 
remnant) and kinematics, yielding distributions of shape and 
kinematic properties (including ellipticity, isophotal shape 
Gil/ a, rotation (v/cr)*, anisotropy, triaxiality, and kinematic 
misalignments) in good agreement with observed cusp or 
disky, rapidly rotating ^ ellipticals. 

The evolution of these sub-systems in re-mergers is not 
readily apparent in the surface brightness profiles of the rem- 
nants - their contribution to the profile is much less than 
the smoothing effects and variations introduced by the re- 
merger (shown in §[4]i. However, their evolution in re-mergers 
is briefly discussed in § [8] and will be studied in detail in 
ICox et all (l2009l) . who show that they are generically de- 
stroyed by major spheroid-spheroid re-mergers. This explains 
the significant effects on the global kinematic properties of re- 
merged ellipticals shown in §[8]- making them rounder, box- 
ier, and more slowly rotating, in agreement with the observed 
properties of massive core or boxy, slowly rotating ellipticals. 
But effects on the velocity field in detail (as a function of ra- 
dius, at higher order in the velocity moments, and in spatial 
distribution rather than just azimuthal average) will be more 
pronounced, and present opportunities not just to test whether 
or not re-mergers are good analogs to boxy/slowly-rotating 
ellipticals, but to distinguish the extent to which substruc- 
ture and kinematic subcomponents contribute to these rota- 
tion/isophotal shape properties, as opposed to e.g. global an- 
gular momentum content or features of the spheroid poten- 
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tial. Kinematic misalignments, kinematically decoupled sub- 
systems, triaxiality, and trends of isotropy w ith radius are also 
more sensitive to orbi tal param eters (see e.g.lCox et al.l2006bl: 
iBovlan-Kolchin etaT 2006; Je sseit et"aLll2007h and may dis- 
tinguish merger histories with preferences for certain orbital 
configurations (possible if e.g. massive galaxies accrete mi- 
nor companions preferentially along filaments). These more 
detailed, "second-order" structural parameters may therefore 
represent a more sensitive probe of the re-merger history (i.e. 
the history of subsequent mergers after the original, spheroid 
forming gas-rich merger), while the "first-order" parameter 
(the light profile) retains a memory of the degree of dissipa- 



tion imprinted in the original gas-rich merger. 
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TABLE 1 
Fits to Core Ellipticals 



Name 


Ref. 


Nphot 
("31 


M* 


My 


(T 

("61 


Re 
\< } 


e 

\°) 


10004/0 


(v/a)* 


% (fit) 


NGC 0584 


2,3 


5 


11.27 


-21.91 


111 


3.39 


0.30 


1.50 


1.55 


2 02+005 

-0 19 


NGC 0720 


2 


4 


11.64 


-22.20 


247 


6.31 


0.39 


0.35 


0.23 




NGC 0741 


2 


1 


12.25 


-23.27 


293 


13.18 


0.15 


- 


- 


3.10 


NGC 0777 


3 


3 


11.65 


-22.92 


348 


7.17 


- 


-0.20 


0.28 


2 '30+0.16 
'^•■'-'-0.16 


NGC 1316 


2 


1 


12.27 


-23.32 


250 


7.24 


0.37 


1.00 


0.91 


2.11 


NGC 1374 


2 


1 


10.71 


-20.57 


207 


2.57 


0.11 


- 


- 


4.37 


NGC 1399 

NGC 1407 


2 

3 


4 

2 


11.56 
11.31 


-22.07 
-22.20 


359 

285 


4.17 

3.32 


0.12 
0.05 


0.10 

-0.20 


0.25 

0.84 


2 £1+0.08 
-0 11 
1 22+0.26 

-Q 26 


NGC 1600 
NGC 1700 


2,3 
2 


4 
4 


11.93 
11.49 


-23.21 
-21.95 


321 

234 


12.96 
4.37 


0.33 
0.29 


-1.20 
0.90 


0.05 
0.80 


1 02+0.26 

-0 31 
4 44+1.04 
^•^ -0 57 


NGC 2832 


2 


2 


12.52 


-23.76 


330 


31.62 


0.30 


-0.30 


0.12 


2 Q^+0.04 
^•'-^-0.04 


NGC 3379 


2,3 


7 


11.03 


-21.14 


221 


2.45 


0.10 


0.20 


0.82 


2 54+0.32 
-0 58 


NGC 3607 


3 


3 


10.88 


-21.22 


248 


1.68 


0.16 


- 


0.92 


1.80+00^ 

2 f,s+0:20 


NGC 3608 


2 


3 


11.02 


-21.12 


195 


3.31 


0.18 


-0.20 


0.44 


NGC 3640 


2,3 


6 


11.17 


-21.53 


180 


3.75 


0.21 


-0.20 


1.48 


2 1 +1.14 
■-^ -0 17 


NGC 3842 


2 


1 


12.19 


-23.18 


316 


17.78 


0.17 


- 


- 


3.09 


NGC 4168 


3 


3 


10.91 


-21.45 


182 


3.94 


0.11 


- 


0.26 


2.55+''-22 


NGC 4261 


1,3 


4 


11.41 


-22.33 


294 


11.09 


0.21 


-1.30 


0.10 


2.40:|| 


NGC 4278 


2,3 


8 


10.70 


-20.52 


259 


2.02 


0.16 




0.71 




NGC 4291 

NGC 4365 


2,3 

1,2,3 


4 

7 


10.71 

1 1 .59 


-20.63 
-22.44 


259 

269 


1.60 

13.83 


0.25 
0.24 


-0.40 

-1.10 


0.52 

0.08 


1 07+1.84 
^•Di-0.44 


NGC 4374 


1,3 


4 


11.10 


-22.43 


287 


8.29 


0.14 


-0.40 


0.09 


2 12+'" 
^•^ -0.34 


NGC 4382 


1,2,3 


4 


11.26 


-22.43 


196 


8.43 


0.19 


0.59 


0.33 


3 89+' " 

-1.35 
2 gj+0.90 


NGC 4406 


1,2,3 


7 


11.36 


-22.66 


250 


19.72 


0.29 


-0.70 


0.18 


NGC 4472 


1,2,3 


7 


11.70 


-23.12 


287 


16.34 


0.16 


-0.30 


0.47 


2 4c+8:S 

-0 95 


NGC 4473 


2,3 


6 


10.91 


-20.82 


178 


2.56 


0.43 


0.90 


0.40 


2 cQ+0.32 
-0 75 
2 70+1.93 

'°-Q 22 

1.06^:1 


NGC 4486 
NGC 4489 


1,2 
1,3 


4 

3 


11.93 
9.82 


-22.95 
-18.56 


360 
49 


31.44 
1.10 


0.12 
0.12 


0.01 
-0.20 


0.11 
1.48 


NGC 4552 
NGC 4589 
NGC 4636 


1,2,3 
2,3 
1,2,3 


7 
6 
8 


11.25 
11.00 
11.38 


-21.46 
-21.45 
-21.96 


261 
215 
191 


7.74 
4.38 
16.28 


0.07 
0.20 
0.18 


0.01 
- 
-0.10 


0.28 
0.57 
0.25 


2 <o+Ll8 

-0 33 

34Qai.« 

2 7 1 +8:'5t 

^•'^-0.89 


NGC 4649 


1,2,3 


7 


11.73 


-22.62 


341 


10.52 


0.17 


-0.50 


0.46 


2 77+0.66 
-1.00 


NGC 4874 


2 


2 


12.37 


-23.49 


290 


61.66 


0.09 


-0.30 


0.22 


2 oa+0.01 

° -8 01 
1 co+0.56 

-0 70 


NGC 4889 


2,3 


4 


12.20 


-23.62 


381 


17.67 


0.35 


-0.25 


0.05 


NGC 5061 


2 


1 


11.53 


-22.01 


194 


3.55 


0.04 


- 


- 


3.55 


NGC 5322 


3 


3 


11.34 


-22.21 


224 


4.68 


- 


-0.90 


0.26 




NGC 5419 


2 


1 


12.30 


-23.37 


315 


15.14 


0.21 


- 


- 


3.78 


NGC 5490 


3 


3 


11.35 


-22.13 


301 


6.31 


0.19 


- 


0.27 




NGC 5557 


2,3 


4 


11.31 


-22.28 


260 


5.94 


0.21 


- 


0.12 


2.78:»«' 
4 87+^'^' 


NGC 5576 
NGC 5813 


2 

2,3 


3 
6 


11.13 
11.20 


-21.31 
-21.68 


188 
238 


3.80 
9.82 


0.30 
0.16 


-0.50 
0.01 


0.22 
0.51 


NGC 5982 


2 


1 


11.50 


-21.97 


250 


5.01 


0.27 






3.02 


NGC 6166 


2 


1 


12.55 


-23.80 


300 


67.61 


0.28 


- 


0.08 


2.73 


NGC 6702 


3 


3 


11.45 


-22.26 


182 


6.39 


0.23 


-0.40 


0.18 


3.11+"-5' 
1 74+8:S 
-0.20 


NGC 7052 


3 


3 


11.48 


-22.46 


275 


9.77 


0.45 


0.01 


0.34 


NGC 7385 


3 


3 


11.95 


-23.53 


259 


10.27 


0.12 


0.01 


0.16 


1 00+1.04 

-0 37 
7 os+0.39 


NGC 7619 


2,3 


4 


1 1 .60 


-22.60 


337 


6.85 


0.24 


0.20 


0.53 


NGC 7785 


2,3 


4 


11.48 


-22.46 


291 


4.56 


0.42 


-1.50 


0.47 


i.75:':« 
2.08!|| 


IC 1459 


2 


3 


11.81 


-22.51 


311 


4.47 


0.26 


- 


0.22 


IC 4296 
IC 4329 


3 
2 


3 
2 


11.75 
12.63 


-23.18 
-23.95 


323 
270 


5.86 
41.69 


0.10 
0.15 


0.00 
- 


0.64 
- 


1 01+0:01 

-Q05 

2 64+0.65 


A0076-M1 


2 


1 


12.07 


-22.96 


- 


5.95 


- 


- 


- 


3.15 


A0119-M1 


2 


2 


12.67 


-24.01 


- 


45.71 


- 


- 


- 


2 30-*-°'°^ 

-0 03 
-"•"'-0.22 


A0168-M1 


2 


2 


12.16 


-23.12 


- 


31.62 


- 


- 


- 


A0193-M1 


2 


1 


12.61 


-23.91 




63.91 








3.74 


A0194-M1 


2 


1 


12.11 


-23.03 


- 


23.84 


- 


- 




2.94 


A0195-M1 


2 


2 


11.89 


-22.65 




14.23 








T 40+0.06 
•'•^^-0.06 


A0260-M1 


2 


1 


12.36 


-23.47 




31.67 








2.64 


A0295-M1 


2 


2 


12.15 


-23.11 




39.81 








3 05+" "' 
■'•"■^-0.01 


A0347-M1 


2 


1 


12.01 


-22.85 




26.47 








2.75 


A0376-M1 


2 


1 


12.43 


-23.60 




47.86 








3.50 


A0397-M1 


2 


1 


12.33 


-23.42 




33.88 








2.72 


A0496-M1 
A0533-M1 


2 
2 


4 
4 


12.82 
11.91 


-24.28 
-22.68 




58.88 
21.38 








4C+O.14 


A0548-M1 


2 


4 


11.95 


-22.75 




20.89 








3.80-1 



ris (sim) 
(12) 



/. (fit) 
(13) 



fsb (sim) 
(14) 



2 77+0.45 
^- ' -0.20 

2 ^Q+1.07 
^■-^ -0.34 
2 41+0.69 

2 no+?:§I 

-Q.S2 
2 48+0.75 

'■^-'-0.35 

o 07+0.44 

2 CO+0.7I 

3 50+'' *^ 

295^1:9^ 

^■^ -0.79 
3 19+1.48 

3 11+'^*^ 

-0.87 

3.78!i« 
2 02+ 1-05 

3.69;"-''^ 

2 go+V^i 

7 70+?:^?) 
^- ''^-0.38 

3 32+'*'^^ 
2 92+A:^5 

7 o-i+l'Sl 

2 57+2-" 
^■■^'-0.35 

3 1 1+0.63 

-0.80 
3 09+' " 

o2+''-91 

2 yi+8:§^ 

3 06+^^^^ 

^■'^°-0.76 
2 79+1.46 
^- ' -0.83 
2 72+'-^^ 

2 91+8:^^ 
^■^'-1.31 

a 02+105 
-0.86 
4.22+"'*'' 

3 10+^-^^ 

3 07+* '2 
^■'-"-0.90 

1.88+2 «2 
2 09+0:33 
^■''^-0.32 

2 75+1.92 
^- '-^-0.79 

3 01+^-33 

2 96+«^^* 
-0.43 

2 03+1.32 
^■^ -0.84 

2 77+0.45 

3 09+^-^^ 

3 ^,+155 
^■■^ -2.83 

1 55+2-25 

2 67+0.64 

^■"'-1.03 

2.53+3-i!« 



-All 



4.67 
2.71 
3.19- 



3.07 

2 91+ii:^^ 

o 10+8:^5 
J.i8_(m 

3 76+0.9 1 

^■'°-1.12 

2 76+2-16 

2 91+^:1')? 

2 09+8:?? 
^■''^-0.36 

3.10:1^ '2 

2 77+8:^^ 

2 77+8:^? 
^- ' -0.20 
2 77+0.45 
^- ' -0.20 

2 U^-^^ 
^•°^-n.4n 



004 
004 
004 
052 



004 
116 
060 



052+" 
n 177+0.043 
"•^ ' '-0.010 

0.026 
n 1 1 0+0.006 

0.075 
0.051 

0.189™ 
n 224+'' 

0.01 1!*!™^ 

0.08 1+j] 
0.124+" 
0.126+^ 
0.108+[j 
0.098+1] 
0.045 

0.081!":""" 

0.105ii 
0.119t!]l? 
0.182!«:»| 

A n67+fl-037 

126+""" 

034+"-008 
"•'"^-0.010 
033+0-011 

091+0-010 

6^^8:8!? 

-S-83^ 
054+0-031 

214+"-02l 
"•^^ -0.022 
136+0-025 

A Q45+0.0i5i 

0.056ii 

0.083!«;M^ 

A Q36+0.000 
U.UJD „ „„„ 

145+*-"" 

"•'^ -0.026 
0.075 

0.086!S;»1| 

0.049 



0.108+[; 
0.088+^ 
0.240+ 



0.044: 



006 
052 



Ml 

-0.001 

0.217 
0.074 

A no2+0-027 

7;8:8V? 

^'•*J^'-0.009 
A 107+"-" 17 
"-"-0.017 

105+"""^ 

"•""-0.006 

A 092+0-006 
-fi.Q26 

A 1 05+0.073 

139+0-003 
-Q.QQ5 
045+0-024 

0.086 



0.025 



+0.001 



0.005!^:»^| 

0.037 

0.056 

060+" ""2 
"•"""-0.002 

0.067 

035+" """ 
"•"^■^-0.000 

0.127 
0.032 

0.079 

030+0-005 

0.034+: 
0.009 



mi 
Mi 



103+"- 191 
"•""-0.066 

096+" '31 
"•"^"-0.061 

076+" "3" 
"•"'"-0.040 

0.076+ 

0.105: 

0.194_Qgj2 

o.io8!«:'| 

0.103! : 2 

0.087!: 
i44+0.1g3 
"•^^^-0.055 

070+" ""'* 

"•"'"-0.025 

190+"'" 

"•^^"-0.092 
229+"- 1 8" 

167+"- '52 
"•'"'-0.077 
154+0.192 



mi 



0.085 
0.077 

0-105:0.068 

230+"-"* 
"•^■'"-0.106 

"•'"^-0.053 
215+"-l™ 

0.098!":"! 

045+0-088 

"•"^^-0.016 

216+"-"'5 
u.zio „ III 

0.083!-' 

0.223!«]i 

190+O-214 

Jio?^8:??i 
"•""-0.Q66 

O76+0-064 
"•"'°-0.043 
103+0-116 
"•""-0.070 
074+0-093 

073+"-"'5 

"•"'^-0.028 

1 94+0-143 
"•'^^-0.103 

o^io3!«:!Jl» 

0.087!!]:!J« 

106+" 121 
"•'"°-0.062 

105+"i"i 
"•'"^-0.068 

0.192 
0.102 

0.118: 

0.089 
0.089 
0.085: 
0.087: 
0.103 
0.046 
0.126: 
0.085: 
0.079 

0.117: 

047+Kii^* 
-9MX 

047+0-058 
"•"^'-0.001 
076+" "*2 

0.074^:| 

0.081« 

060+"-"^' 
"•"""-0.014 

0.047!":"! 

0.087!: 

0.047!«:i 

0.084!S:«^ 

0.065+"-"?5 



mi 
m\ 
m^ 

J:?S 



mi 

Mil 

mt 



Mi 



0.085 
0.077 



mi 



44 



Hopkins et al. 

TABLE 1 — Continued 



Name 


Ref. 


Nphot 


M* 


Mv 


a 


Re 




(v/aY 


n.! (fit) 


n, (sim) 


/. (fit) 


fsb (sim) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) (9) 


(10) 


(11) 


(12) 


(13) 


(14) 



A0634-M1 


2 




11 


92 


-22 


70 




16.60 








3.24 


A0779-M1 


2 




12 


54 


-23 


78 




36.23 








3.01 


A0999-M1 


2 




11 


78 


-22 


45 




13.80 








2.74 


A1016-M1 


2 




11 


67 


-22 


26 




22.41 








4.82 


A1142-M1 


2 




12 


01 


-22 


85 




24.84 








4.99 


A1177-M1 


2 




12 


41 


-23 


56 




66.74 








4.31 


A1314-M1 


2 




12 


27 


-23 


31 




96.02 








4.34 


A1367-M1 


2 




12 


12 


-23 


05 




14.86 








2.88 


A1631-M1 


2 




12 


29 


-23 


34 




32.36 








4.28 


A1656-M1 


2 




12 


55 


-23 


81 




22.55 








2.20 


AZU4U-M 1 


z 




1 1 
12 




—li 


40 












1.04 


A2052-M1 


2 




12 


11 


-23 


04 


250 


189.2 


0.24 






3.48 


A2147-M1 


2 




12 


18 
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Compiled and fitted parameters for the confirmed core ellipticals in our observed samples. Columns show: (1) Object name. (2) Source for surface brightness 
profiles, where 1 MKormendv et alj 420081) , 2 = (Lauer et alj 42007 al) , 3 = (Bender et alj 09881) . (3) Total number of different surface brightness profiles in our 
combined samples for the given object. (4) Stellar mass [logM* /Mq]. (5) V-band absolute magnitude. (6) Velocity dispersion [km/s]. (7) Effective (half-light) 
radius of the total light profile [kpc]. (8) Ellipticity. (9) Boxy/diskyness. (10) Rotation. (11) Outer Sersic index Us of the two-component best-fit profile. 
Where multiple profiles ai^e available for the same object, we show the median and minimum/maximum range of fitted values. (12) Range of outer Sersic 
indices fit in the same manner to the best-fit simulations, at f fa 1 — 3 Gyr after the merger when the system has relaxed. (13) Fraction of light in the inner or 
extra light component of the fits. Where multiple profiles are available for the same object, we show the median and minimum/maximum range of fitted val- 
ues. (14) Fraction of light from stars produced in the central, merger-induced starburst in the best-fit simulations (± the approximate interquartile range allowed). 
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Fig. 33. — As Figure [Tol but in log-log space. Surfa ce b rightness profiles - decomposed into dissipationless and dissipational/extra light components — are 
shown for core ellipticals in and around the Virgo clustefiSl. Open circles show the observations, from lKormendv et alj 120081) . These are the highest-mass core 
ellipticals in Virgo (~ 6 — 8M»). Upper: Observed V-band surface brightness profile with our two component best-fit model (solid, dashed, and dotted lines 
show the total, inner/extra light component, and outer/pre-starburst component). The best-fit outer Sersic index, extra light fraction, and rms residuals about 
the fit ar'e shown. Lower: Colored lines show the corresponding surface brightness profiles from the three simulations in our library which con'espond most 
closely to the observed system. Dashed line shows the profile of the starburst light in the best-matching simulation. The range of outer Sersic indices in the 
simulations (i.e. across sightlines for these obj ects ) and range of starburst mass fractions which match the observed profile are shown, with the rms residuals 

of the observations about the best-fit simulatiotffH Bottom: Observed disky/boxy-ness (04) and ellipticity profiles, with th e media n (solid) and 25 — 75% range 
(shaded) corresponding profile from the best-fitting simulations above. Note that these ar'e not fitted for in any sense. Figures [34|37| show the other core ellipticals 
in the sample, ranked from most to least massive. 

APPENDIX 

FITS TO THE SAMPLE OF KORMENDY ET. AL. 2008 



In Figures l33]l37l we reproduce Figures [T0lfT4l but with profiles shown in log-log projection as opposed to r'/^ projection. 
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Fig. 34. — The next most massive core ellipticals (~ 5 — 6M»). (As Figure [TTI but in lo; 
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Fig. 35. — Similar massive core ellipticals (~ 4 — 5M,). (As Figure [T2] but in log-log space.) 
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The next most massive core ellipticals (~ 2 — 3M»). (As Figure [73] but in log-log space.) 




